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Abstract 



This thesis describes some of the recent (and some less recent) developments 
in calculational techniques for scattering amplitudes in quantum field theory. 
The focus is on on-shell recursion relations in complex momenta and on the use 
of unitarity methods for loop calculations. In particular, on-shell recursion is 
related to the MHV rules for computing tree-level gauge amplitudes and used 
to extend the MHV rules to graviton scattering. Combinations of unitarity cut 
techniques and recursion are used to argue for the "No- Triangle Hypothesis" in 
Af = 8 supergravity which is related to its UV behaviour. Finally, combinations of 
unitarity and recursion are used to demonstrate the full calculation of a one-loop 
amplitude involving a Higgs particle and four gluons in the limit of large top mass. 
The present version is edited to incorporate some of the comments and suggestions 
of the evaluation committee, but has not been updated for developments in the 
meantime. 
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Chapter 1 

Introduction 



1.1 Scattering Amplitudes in Theory and Experiment 

Ever since Physics arose as an empirical science, one of its main objectives has 
been the identification and understanding of the fundamental principles and con- 
stituents of Nature. Over time, fundamental principles have changed or been 
extended, and constituents regarded as fundamental have been found to contain 
even more fundamental parts. For the last 100 years, this hunt for the fundamen- 
tal has been tied to the principles of quantum mechanics which governs Nature 
at the fundamental level. 

Particle physics as we know it today studies objects which are so small and 
intangible that they can rarely be held in one place for study, either because they 
are massless and move at the speed of light, or because they simply decay too 
fast. At the same time, the measurement apparatus cannot possibly be small 
enough to explore a fundamental particle without being a particle itself. Thus, 
the only viable method of exploration is to study the interactions between pairs 
of particles that collide by identifying the consequences of this collision: What 
comes out where, and how fast. 

The physical quantity that describes this is the scattering cross section. It 
specifies the area within which one particle must hit the other for a specific 
process to take place, or rather, the area where the particles do collide times the 
quantum mechanical probability (density) that the given process will occur if the 
particles collide. Such quantum mechanical probability densities always arise as 
the absolute value squared of a quantum mechanical amplitude, which is in this 
case called a scattering amplitude. 

A scattering amplitude can receive contributions from different processes 
which are indistinguishable in the quantum sense, that is, processes which look 
the same to the measurement apparatus. To form a scattering cross section, all 
these processes have to be summed before squaring. Thus, if we wish to say 
something meaningful about a theory with certain particles but not all of those 
we know to exist, we should not start squaring the amplitude as we would be 
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missing interference terms from the other processes. This point is made to argue 
that the scattering amphtude is the furthest one can calculate in a theory without 
taking into account other physical processes or experimental conditions. In that 
sense, the scattering amplitude is the most proper way to describe a result in any 
reductionist approach to particle physics. 

Scattering amplitudes arc^ normally calculated by perturbation theory as power 
series in coupling constants. The lowest order results are easily accessible with 
today's methods, but already the next order is an analytical and numerical chal- 
lenge. This is unfortunate because, as we shall see, more precise (approximate or 
exact) results are in high demand. 

1.1.1 Background Processes at the LHC 

For the next many years, particle physics will be dominated by the experimental 
programme at CERN, where the Large Hadron Collider (LHC) is being built. 
This machine will collide (for the most part) protons on protons with a center-of- 
mass energy of 14 TeV. This is hoped to produce an array of novel particles, such 
as supersymmetric partners of all the known particles, and expected to produce 
decisive statistical evidence for the existence of the Higgs boson. 

These interesting results will, however, be hiding behind massive amounts of 
particle reactions that we already know and love from earlier investigations of 
QCD and electroweak interactions. A good knowledge of the production rates 
of these well-known events is essential because they contaminate the interesting 
events and need to be subtracted. Knowing only these background processes to 
the lowest order in the coupling constants is likely to introduce large theoretical 
uncertainties, some of them prohibitively large. 



process 

{Ve{Z,W,^}) 


relevant for 


1. pp — >■ yyjet 

2. pp tt bb 

3. pp ^ tt + 2 jets 

4. pp^VVbb 

5. pp^ -1-2 jets 

6. pp ^ V + 3 jets 

7. pp^VVV 


tiH, new physics 

tiH 

tiH 

VBF^ H VV, tiH, new physics 
VBF^ H ^VV 
various new physics signatures 
SUSY trilepton 



Table 1.1: The LHC "priority" wishlist, extracted from [65]. 



Theorists and phenomenologists working on the LHC have issued a prioritized 
wishlist of next-to-leading order calculations they want done before the LHC 
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turns on [65] (scheduled for Spring 2008 at the time of writing) included as table 
11.11 At the time of writing, only the first and the last seem to be close to 
solved [70,75,106]. In principle, it is known how to calculate all these, but the 
conventional methods seem to have hit a roadblock; the calculational complexity 
of the required calculations is rising faster than the available resources, be they 
human or electronic. Decisive progress in this area must come from both improved 
methods of calculation and a strengthened effort. 

1.1.2 Approaching Quantum Gravity? 

Apart from the prosaic undertaking of extracting useful results from a multi- 
billion Euro experiment, scattering amplitudes also play a role in the theoretical 
understanding of quantum field theories. Amplitudes express and reveal both the 
internal symmetries of a theory and its internal inconsistencies. 

One inconsistency plays a special role in the problems defining a quantum 
theory of gravity, namely that of non-renormalizability. In contrast to other 
quantum field theories where infinities that pop up in intermediate results can be 
arranged to vanish in physical quantities, gravity theories normally require more 
and more new terms in their defining equations to counter the infinities which 
necessarily appear. 

However, it so happens that the maximally supersymmetric {M = 8) quan- 
tum field theory of gravity does not seem to have such infinities. Calculation 
of the necessary scattering amplitudes in this theory (and others) of gravity are 
unfortunately extremely cumbersome and evidence is only accumulating slowly 
that this may indeed be the case. 

This understanding is driven by new methods of calculating amplitudes. 
Thus, it may just be that such new methods can give a glimpse of something 
which, without being a true description of all known particles and interactions, 
can be called a consistent quantum theory of gravity. 

1.2 Scattering Amplitudes in the Complex Plane 

Scattering amplitudes are usually calculated perturbatively using Feynman rules. 
Feynman rules are derived directly from an action principle, are understood by 
all particle physicists alike, and have well-studied mathematical properties. A 
calculation done with Feynman rules is rarely called in question, save for the 
pointing out of simple algebraic errors. This makes Feynman rules the preferred 
choice of almost everyone in the community, and with good reason. 

Unfortunately, there are some issues with Feynman rules that make them un- 
fit in several situations. The first is the issue of complexity: Since the rules work 
by summing over contributions from all graphs of certain kinds, the complexity of 
the calculation is factorial in the number of participating particles. When going 
beyond leading order in coupling constants, there are integrals whose complex- 
ities are also factorial in the number of particles. This leads to the roadblock 



7 



mentioned earlier. The second issue is that of unphysical singularities: The re- 
sults of the calculation are most often quoted in a form where some singularities 
have to cancel between terms. When evaluating such expressions numerically, 
roundoff error may introduce apparent effects which are not really there. 

A contributing factor to these problems is that Feynman rules throw out in- 
formation in exchange for more mechanical calculations. The prime example of 
this is non-abelian gauge theories, where gauge invariance enforces a particular 
structure of the three and four-point interactions. In the actual calculations, how- 
ever, no reference is made to the gauge invariance, which might as well not have 
been there. In the end, the constrained form of the Feynman rules ensure that 
the amplitude is gauge invariant, but this has been obscured by the intermediate 
calculation. This is just one of several pieces of information which may be used to 
circumvent Feynman rules and obtain amplitudes faster, in more compact forms, 
and with fewer unphysical singularities. 

1.2.1 The Unitary and Analytic S'-Matrix 

The focus of this thesis will be the use of methods of complex analysis and meth- 
ods derived from such considerations. Although most of the methods described 
here have been developed since 2004, the idea itself is by no means new. In some 
sense, it can be traced back to the so-called ^-matrix programme of the 60's. 

The S'-matrix is the evolution operator and is conventionally decomposed into 
a unit piece (which describes no scattering) and into the rest (which describes 
the actual scattering). 



The scattering operator iT sandwiched between a state of particles coming in 
from the infinite past and a state of particles coming out at infinite future, gives 
the scattering amplitude (also known as an S-matrix element if the unit matrix 
is included). Requiring that the S'-matrix be unitary imposes a relation for the 
T operator commonly known as the optical theorem. 



If we insert a complete set of states between T and on the right side, we get 
the formal diagrammatic expression 



S=l + iT. 




2imr = rr^ 





(1.3) 
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where the sum is over ah possible insertions of internal states and the integral 
is over their on-shell momenta. In perturbation theory, the left side could cor- 
respond to a one-loop amplitude and the right could correspond to a product of 
two tree amplitudes. 

If we also view the S'-matrix as an analytic function of the kinematical vari- 
ables, we can gain more information. In general, a one-loop amplitude as a 
complex function of one of the kinematic variables will have a branch cut at real 
positive values of the variable. The optical theorem can then be twisted to give 
the discontinuity across the branch cut in terms of the insertions of all possible 
states between two tree amplitudes, and the hope would be to reconstruct the 
one-loop amplitude from its branch cut discontinuities. The really interesting 
thing would then be if this information could somehow be used to bootstrap 
amplitudes to all orders in perturbation theory. 

This was the goal of the S'-matrix programme: to produce a theory for the 
strong interactions (and, in the long run, all of particle physics) from unitarity and 
analyticity constraints on the S'-matrix. That programme failed, partly because 
a host of other assumptions turned out to be needed for the full reconstruction of 
the S'-matrix, partly because QCD, a Lagrangian quantum field theory, became 
established as the correct theory of the strong interactions. A review of the 
methods and motivations of the S-matrix programme can be found in [83]. 

In the 90's, it was realized that some of the tools of S-matrix theory were 
more powerful than public opinion in the community had assumed. The power, 
however, came from the input of facts known from Feynman rules and dimensional 
regularization, two concepts which were absent in the original S-matrix theory. 
Bern, Dixon, Dunbar, Kosower and others used the known structure of certain 
one-loop amplitudes to write them in a form where the discontinuities across 
branch cuts could be expressed easily so that a (not even complete) computation 
of the right side of (II. 3p would almost completely determine the amplitude [26, 
33-35] . We return to a detailed describtion of these methods in chapter [5l 

Although this so-called "unitarity cut technique" produced results which were 
unthinkable to compute with Feynman graph techniques, and attempted to use 
all known information as efficiently as possible, it did lack in generality only being 
useful for some particular cases (some of which were nevertheless relevant to ex- 
periment). A central point in these methods is that only on-shell tree amplitudes 
are used as input for the method, and this hints that we may, after all, be able 
to bootstrap our way through the perturbative expansion. 

This is not as far fetched as it may sound, because there are already results 
that indicate that this is feasible in principle. The Feynman Tree Theorem [87, 
100] gives exactly such a construction. It works by using retarded propagators 

in all propagators in a loop, which set the amplitude to zero because of time 
ordering. By using that the retarded propagator is the Feynman propagator plus 
a delta function 

Ar{P) = ApiP) - 2TTi5{P^ - m^)e{-Po) (1.4) 



9 



the expression can be rewritten as a sum of ordinary Feynman diagrams where 
some number of propagators are Feynman and others are set on-sheU by the 
delta function. This construction shows directly that the all-orders perturbative 
S'-matrix can be written as some complicated convolution of the tree-level S- 
matrix. The Feynman Tree Theorem is discussed in more detail in [55]. 

1.2.2 Generalizing from Loops to Trees 

As we saw above, the use of complex momentum space methods is not a new 
invention, and at one-loop level, it has been used for some time as a tool for the 
most efficient evaluation of scattering amplitudes. Most tree-level calculations 
were still being conducted in real Minkowski space using Feynman diagrams (or 
related methods) at the turn of the century. The change in affairs would, however, 
come from a completely different direction. 

It had been known since the mid-80's that a certain class of tree amplitudes 
were considerably simpler than a Feynman diagram calculation would suggest 
when using the spinor helicity formalism, to which we will turn in chapter [21 
These were the so-called Maximally Helicity Violating (MHV) amplitudes where 
(when all gluons are considered outgoing) two gluons have negative helicity and 
any number have positive helicity. Building on some critical earlier insight by 
Nair [115], Witten [133] was able, in 2003, to describe gluon tree-amplitudes 
from the topological B-model string theory in twistor space, a type of theory 
not commonly associated to the nuts-and-bolts of scattering calculations. Soon 
after it was realized that in the same sense that MHV amplitdes are lines in 
twistor space, tree-amplitudes with more negative helicity gluons were collections 
of intersecting lines in twistor space. From this, it was only a small leap to 
conjecture that there existed a formalism for calculating gluon amplitudes whose 
vertices were MHV amplitudes. 

That this was indeed the case was shown by Cachazo, Svrcek and Witten in 
the beginning of 2004 [69]. The construction is commonly known as either MHV 
rules or CSW rules. This thesis will use the first, knowing that it may cause a bit 
of confusion (Feynman rules generate Feynman amplitudes; MHV rules do not 
just generate MHV amplitudes). A description of MHV rules and some extensions 
can be found in chapter [3l A good deal of calculations followed which verified 
the MHV rules [91,102,135,136]. Brandhuber, Spence and Travaglini were then 
able to use them for one-loop calculations [57] which was slightly surprising at 
the time. 

All of these events were directly related to the twistor picture of gauge ampli- 
tudes, which in itself did not involve complex momenta, but to some extent it did 
involve momenta in 2 + 2 dimensions. In that metric signature, Britto, Cachazo 
and Feng [59] realized that the unitarity cut technique described above can be 
extended to insert four on-shell internal states rather than just two, basically 
because the on-shell demand is easier to satisfy with two time directions. When 
there are four on-shell constraints in a loop integral it becomes trivial because 
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the four coordinates of loop momentum are fixed. This so-called "Quadruple Cut 
Technique" allowed e.g. the conversion of the calculation of one- loop amplitudes 
in maximally supersymmetric {M = 4) Yang-Mills to a purely algebraic one in 
terms of products of tree amplitudes. It was also realized that the constraint 
of being in 2+2 dimensions could effectively be removed by assuming to be in 4 
complex dimensions. We return to this technology in chapter [5l 

Since there are relations between tree and one-loop amplitudes in = 4 
super- Yang-Mills [31,71,92,105], this allowed a quartic recursion relation for on- 
shell tree- level Yang-Mills amplitudes [128], but with complex momenta. Britto, 
Cachazo and Feng [60] were then able to rewrite this as a quadratic recursion rela- 
tion, but they also showed together with Witten [61] that it resulted from taking 
the external momenta seriously as complex and using simple complex analysis. 
Moreover, they showed that these recursion relations could be iterated such that 
only complex three-point amplitudes contributed; in other words, the four-point 
vertex in Yang-Mills can be thought of purely as an artifact of gauge invari- 
ance necessary for Feynman rule calculations. These results will be presented in 
chapter [3l 

The remarkable developments described here have presented a completely new 
view on the calculation of scattering amplitudes, and tie together concepts across 
loop order in ways which are still not properly understood. The consequences of 
treating scattering amplitude calculations as complex analysis is, at its essence, 
the topic of this thesis. 

1.2.3 Full Amplitudes from Trees 

The concept of on-shell recursion relations derived from complex-momentum 
methods have numerous applications. Most obvious are concrete tree-level cal- 
culations which come out in a form which is more compact than the form from 
Feynman rules. Another tree-level application is to show how the MHV rules 
mentioned above can be viewed from recursion as a particular reorganization of 
the Feynman rule calculation. The connection between MHV rules and on-shell 
recursion will be explored in detail in chapter [H 

Another use of on-shell recursion is for the calculation of one-loop amplitudes. 
Parts of one-loop amplitudes can be computed using the unitarity cut technique, 
and exactly the parts that are missed in that technique are rational ones which 
have similarities to tree amplitudes. This means that they are also computable 
by recursion relations [38-40] although in a slightly more involved form. Such 
a determination of a one-loop amplitude from unitarity cuts and recursion is 
presented in chapter [51 

Numerous developments of unitarity combined with recursion have accumu- 
lated over the last years and have improved one-loop calculation. Interestingly, 
some results extend to a higher number of loops, although they are primarily 
unitarity based. Founded on earlier insight by Anastasiou, Bern, Dixon and 
Kosower [2], Bern, Dixon and Smirnov [41] have been able to propose a closed 
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form result for the all-foops AA = 4 four and five-point amplitudes. Similarly, 
there have been calculations up to three loops in M = 8 supergravity which 
lends credibility to the notion that it is perturbatively UV finite, and hence a 
(comparably) well-defined theory of gravity outside string theory. The input into 
the calculations are tree amplitudes and the results are consequences of their 
structure. These multi-loop results continue to attract attention from outside 
the circle of specialists. 

1.3 About This Thesis 

1.3.1 Omissions 

The purpose of this thesis is to review some of the central complex momentum 
methods developed in recent years and to convey the results of the authors own 
works in this area. As the reader will notice, the review part already takes up a 
large part of the thesis, but still the topics that are relevant and interesting in 
this context are so numerous and complex that many of these things have had to 
be omitted. 

Twistors are omitted completely, a choice which could be percieved as radical 
when the field has such close connections to twistor string theory. In the approach 
taken here, however, this connection is more historical since complex momentum 
methods has managed to outrun the twistor methods in many cases, in particular 
those cases chosen here. Twistor string theory has been well reviewed by Cachazo 
and Svrcek [67] to which the reader is referred. Gauge theories in twistor spaces 
have been developed by e.g. [52]. 

Another area which is not covered is that dealing with MHV rules as derived 
from ordinary Yang-Mills in light-cone gauge. In this approach [93,113], the light- 
cone Lagrangian is rewritten in a form where the MHV rules appear naturally, 
and this permits extensions to loop calculations. Some important works in this 
area are [56,86]. The viewpoint taken on MHV rules here is that MHV rules are 
a special case of recursion relations, a viewpoint which is admittedly too narrow. 

There have also been refinements of many of the unitarity and recursions 
mentioned here which would lead us too far if they were to be explained in 
detail. These are mentioned in the relevant connections. Within more traditional 
calculational approaches, there have also been recent advances which will not be 
covered. 

1.3.2 Outline 

The thesis is roughly divided in two parts, having to do with tree-level and 
one-loop level, and is not in historical order. It starts out by explaining some 
basic results in Yang-Mills theory, some notation, and the concept of gravity as a 
quantum field theory in chapter [2j Armed with those methods, we can take a look 
at the two important tree-level developments, namely the MHV rules and on-shell 
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recursion in chapter[3l Based primarily on the articles [126] by the author and [49] 
together with Bjerrum-Bohr, Dunbar, Ita and Perkins, chapter [J] describes how 
MHV rules can appear as a special case of recursion relations, and how that can be 
used to derive MHV rules for a gravitational theory. This highlights some of the 
similarities between gravity and Yang-Mills that are obscure in their Lagrangian 
formulation. 

The purpose of chapter [5] is to introduce the loop level methods that will guide 
us through the rest of the thesis, as well as some related ones, and the remaining 
three chapters present some more concrete one-loop calculations which serve dif- 
ferent purposes. The calculations of chapter [6] serve to explain how amplitudes in 

= 4 super- Yang-Mills are computed using a variety of techniques. The part of 
the calculation which regards amplitudes with other external field content than 
gluons is based on [127] together with Bidder and Perkins. Chapter [7] contains 
calculations which elucidate the one-loop structure of maximal {M = 8) super- 
gravity, which is necessary for the understanding of the all-loop structure. As 
such, it serves to support the notion that M = 8 supergravity is UV finite, and is 
based on [50] together with Bjerrum-Bohr, Dunbar, Ita and Perkins. The thesis 
ends on a more prosaic note with the calculation in chapter [8] of the amplitude 
for a Higgs going to four gluons of particular helicities in the limit where the top 
mass is large. It is based on work together with Badger and Glover [10] , but the 
treatment is slightly different from the article. 

1.4 Conclusion and Outlook 

Rather than ending the thesis with a conclusion, we will sum up the main points 
that the thesis will support, and give an outlook for the future. 

1.4.1 A Competitor in the Calculation Race 

The combined methods of unitarity and recursion have provided some extremely 
effective means of calculating one-loop scattering amplitudes. More specifically, 
they have allowed for the calculation of hitherto unknown six-parton QCD am- 
plitudes, and hold the promise of making the progression to seven-parton QCD 
amplitudes smoother. As demonstrated in chapter [HI they have also allowed for 
calculations with fewer external particles to be sufficiently simplified that com- 
pact algebraic expressions for amplitudes can be generated without the use of 
computer algebra. Thus, it has opened up a new frontier for the specialists, and 
hopefully provided tools for the many smaller groups of physicists with simpler 
calculations in mind. New twists on these methods seem to keep coming in. 

This is not to say that these methods are the only game in town. A purely 
Lagrangian approach is still favoured by the broad community, and there have 
been several refinements of conventional technology, both semi-numerical and 
algebraic, which are reaching a comparable level. One can only hope that a 
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combined effort of all the methods on the market can have the results in place 
when they are needed in the near future. 

1.4.2 Understanding Quantum Field Theories 

The use of complex momentum methods have shown us that perturbative quan- 
tum field theory is still far from understood. The investigation of scattering 
amplitudes as complex functions gives several puzzles. At tree level, the simplest 
version of on-shell recursion described in chapter [3] is known to work for both 
Yang-Mills theories and gravity, but more complicated versions, such as those 
used in chapter HI are not proven but seem extremely likely to hold. Superfi- 
cially, the existence of on-shell recursion seems related to the UV behaviour of 
tree amplitudes, a relation which is not clear either. 

Equally, if not more, interesting are the relations between loop orders. Using 
unitarity, statements about loop amplitudes seem to be reformulable as state- 
ments about tree amplitudes. In chapter [71 we will see how the existence of 
on-shell recursion for tree gravity amplitudes comes close to proving that the 
one-loop structure of = 8 supergravity is the same as = 4 super- Yang- 
Mills. Putting such relations on a firm footing would contribute greatly to both 
our fundamental understanding of perturbative quantum field theory, as well as 
our ability to do calculations. 

Together with the understanding of the full perturbative expansion of Af = 4 
Yang-Mills and the prospects for a UV finite theory of gravity, this makes the 
topic of this thesis most exciting, an excitement that I hope that you, dear reader, 
will share. 

Enjoy! 
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Chapter 2 

Preliminaries 



2.1 Colour Ordering 

One of the more annoying aspects of doing gauge theory calculations is the pro- 
liferation of indices, both for colour and for space-time. A strategy to avoid this 
is to specify quantum numbers (such as colour) as early as possible in the calcu- 
lation. This can provide expressions with less index mess, compensated by more 
expressions to compute. The two chief methods for doing this are colour ordering 
as described in this section and spinor helicity notation described in the next. 
The contents of sections 12.11 to 12.31 are also reviewed in [76] . 

The strategy inherent in colour ordering is to identify the possible colour 
structures we can obtain when we have finished an on-shell calculation. If we 
consider adjoint particles only, the colour structures are products of structure 
constants, /"'"^, which can be written in terms of traces of Nq x Nc fundamental 
colour matrices, 

/«^^ = -^^Tiieth" - m% (2.1) 

where we normalise the matrices as 

[^^^*] = iVir^H^, Tv[tH^] = 6"^. (2.2) 
By using contraction identities such as (for SU{Nc)) 

(t^iit^i = 6iSi - ^6161, (2.3) 

we can write all colour structures in terms of traces of fundamental colour ma- 
trices. Note that when we are contracting structure constants only, the last term 
of (j2.3p will always drop out. Thus, we can write the completeness relation as 

Tr[t''M]Tr[t''N] = Tr[MN] (2.4) 

where M and N are arbitrary matrix strings. 
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At the end of the calculation, (the tree-level part of) an amplitude can be 
written as a linear combination of the distinct trace structures, 

An= Yl Ti[t''-H'''^2...e^njj(^(ai,a2,...,an), (2.5) 

{a,}eS„{l,...,n)/Z„ 

the simplest examples being 

^3 = Tr[t^H''H''^]A3{l,2,3) +Tr[t''H''H''^]A3{l,3,2) (2.6) 

and 

Ai = TT[t''H''H''H''^]A4{l, 2, 3, 4) + Tv[t''H''H''H^^]A4{l, 2, 4, 3) 
+TT[t''H''H''H'''']A4{l, 3, 2, 4) + TT:[t''H''H''H''^]A4{l, 3, 4, 2) 
+TT[t''H°'H''H''^]A4{lA,'^,^) +Tr[t°'H''H''H''-']A4{l,4:,3,2). (2.7) 

The coefficients of the trace structures are called colour ordered amplitudes, be- 
cause they correspond to a particular ordering of the colour matrices. Instead of 
calculating the amplitude, using (|2.1|) to (|2.4|) . and deducing the colour ordered 
amplitudes we can calculate them directly by using colour ordered Feynman rules. 
A diagram then contributes to a colour ordered amplitude if the cyclic ordering of 
the external particles corresponds to that of the wanted colour ordered amplitude. 

If we are dealing with particles in the fundamental representation, the strings 
of matrices will not be traced, because the colour indices of the fundamental 
particles are free. It is still possible to define a colour ordering, e.g., 

A,{ql,q^j,3-^4-') = ^ t'^lj^4(gi, 3, 4, ga) + [t"*t'^^];;^4(9i, 4, 3, 52). (2.8) 

If we work in an SU{Nc) gauge theory, we now have to take the last term of (|2.3|) 
into account. This corresponds to subtracting at each vertex the contribution of 
a U{1) gauge boson. 

At loop level, we will also encounter multi-trace terms which we need to 
compute. At one-loop level, however, it can be shown that all two-trace terms 
are permutations of the ("planar") single-trace terms [33]. 

2.2 Spinor Helicity Notation 

After trading the colour information for more diagrams, we can now trade the 
use of polarization vectors for yet more diagrams. This consists of writing the 
amplitudes completely (or rather, to the extent possible) in terms of Weyl spinors 
of massless particles. We will find that in this notation we can write down explicit 
polarization vectors corresponding to positive and negative helicity gluons. 

We start by deriving some basic facts about Weyl spinors. The Weyl spinors 
corresponding to a massless momentum can be found by solving the massless 
Dirac equation 

P/.l^uip) = Mp) = (2.9) 
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where u{p) is a four component vector and is a vector of 4 x 4 matrices obeying 
{7^7-} = 25^-/. (2.10) 
We use "Peskin and Schroeder" conventions [122] where 

Switching now to Weyl-van der Waerden notation the Dirac spinor can be written 
as two two-component spinors 

«(p)=(;-:s|^j) (2.12) 

whose indices are raised and lowered with eap, e^^, etc., using southwest-northeast 
contraction. We define antisymmetric inner products between spinors, 

(Ar/) = eaf^X'^v'' = e'^'XaVp = A"r/„, (2.13) 
[Afy] = e^^'X'^f = el^^hfj^ = X^t , (2-14) 

thus we are working with the "QCD" sign convention rather than the "string" 
sign convention. In this notation (12. 9|) becomes 

p-ap^u'^{p)=p-a^''ua,{p) = 0, (2.15) 

which can further be elucidated by noting that p ■ o""" can be decomposed into a 
product of two two-component spinors because det(p ■ a) = p^ = 0: 

p.fj"" = A"(p)A"(p), (2.16) 

which is solved by 

X =V\Po+P3\, X =sign(po) (2.17) 

V\P0+P3\ 

Ai = sign(po)-7=^=, A2 = -V\Po+P3\- (2.18) 
VlPo+Psl 

r = sign(po)(A°)* (2.19) 

From the above properties of the two-component spinors it follows immediately 
that the Dirac spinors satisfying the Dirac equation are spanned by 

ip) = ( ^i^^ )^\P), u~ ip) = ( ) = |p] , (2.20) 
whose conjugate spinors are 



n-ip) = ( Q^^^ J ^ {p\, u-^ip) = (^s^.^^^ ) ^ [p\. (2.21) 
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The definitions of \p), \p] etc., ensure that {pq) = {\{p)\{q)) and \pq\ = \\{p)\{q)\ 
thereby helping to simphfy the algebra significantly. In terms of momenta, the 
antisymmetric inner products are 



ijpq) = -sign(po) 



\pq] = sign (go) 



90 +93 



^•0+^3 



(pi+ip2) + sign(go) 



Po +P3 



90 + 93 



(gi+ig2),(2.22) 



90 + 93 



bo +P3 



{pi - ip2) - sign(po 



Po + P3 



90 + 93 



which obey 

{pq)[qp] = {p + qf 



2p ■ q, [qp] = sign(po9o)(P9)*- 



(gi-iga), (2.23) 



(2.24) 



Wc will often encounter these spinors strung together using gamma matrices. 
The following identities can be shown by using standard Dirac algebra: 



(^17^9] = [qb^P), 
{q\r\q] = 29^ 
{p\l^\Q]{rhM = -2{pr)[qs] 

{ph>^r\Q) = -{q\ri''\p), 

{pq){r \ + {qr){p\ + {rp){q\ = 



(Fierz rearrangement), 



(Schouten Identity). 



(2.25) 

(2.26) 
(2.27) 
(2.28) 
(2.29) 



The last two have equivalent forms with {■] — *■ [-1. We will suppress Feynman 

slashes in general unless possibilities of misunderstandings arise. In fact, we 
will quite often write momenta as their slashed versions rather than their vector 
versions, as the whole point of spinor helicity notation is to get rid of explicit 
vectors. When we write massless momenta in matrix form we should really write 



P = \p)\p\ + \p]{p\, 



(2.30) 



but since only one of them avoids being projected out in spinor strings such as 
(g|p|r] or (9|p|r|s), we will allow ourselves the sloppiness of writing 



P=\p)\p\ or p=\p]{p\. 



(2.31) 



We now know how to treat expressions with Weyl spinors and momenta. 
The last objects which may occur in a colour ordered scattering amplitude are 
polarization vectors. We use the two explicit realizations 



V2{qp) ' 



^nip) 



V2\pq] ■ 



(2.32) 
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In both of these, q is some massless reference momentum. Using the rules above it 
is quite straightforward to verify the standard properties of polarization vectors, 



e±(p)-p = 0, e+{p) = {e-{p))\ (2.33) 

e+(p)-e-(p) = -1, e+(p)-e-(p) =0. (2.34) 

We can argue that q is an arbitrary (apart from p ■ q ^ ^) gauge vector by 
considering the difference beween two different choices of q: 

e^{p;q')-e;{p;q) (2.35) 

_ {Q'hM {qIiM 



V2{q'p) V2{qp) 



(2.36) 
(2.37) 



V2{q'p){qp) 
V2{qq' 



V2{q'p){qp) 

p" (2.38) 



{q'p){qp) 



Pf.. (2.39) 



Thus, the two choices give the same results due to the Ward identity. Finally, 
we note that the polarization sum is that of a lightlike axial gauge with gauge 
vector q: 

+ S4 = + (2.40) 

This particular gauge is ghost-free, which will extremely helpful in later chapters 
when we deal with loops. 

Recently, a Mathematica package for manipulations with spinors in this con- 
text called S@M has been published [109] , but using slightly different conventions. 



2.3 Tree-Level Structure of Yang— Mills Theory 

This section collects some facts of Yang-Mills theory that we will need in the 
future, and is not intended as a review of the subject. We will return to the 
one-loop structure of Yang-Mills amplitudes in section 15.11 

2.3.1 MHV Amplitudes 

It has already been mentioned that Yang-Mills scattering amplitudes take on 
very simple forms in spinor helicity notation. The most striking example of this 
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was conjectured by Parke and Taylor [120] and proven by Berends and Giele [16], 
and states that 



^(l+,...,n+) =0, 
^(l+,...,i-,...,n+) =0, 

•^(^"•••'-"'-'^''•••'-")- (12)(i!)!'.(nl) ' ^'-''^ 

where . . . denotes any number of positive hehcity gluons. The last of these is 
called the "Maximally Helicity Violating" amplitude, or MHV amplitude for 
short. There is a simiar expression for the amplitude with two positive and 
the rest negative helicity obtained by complex conjugation, 

"'' = '-'>"' lizipa^.i,.!! ' <^^^2> 

called the "googly MHV amplitude''^. These formulas determine all gluon ampli- 
tudes up to and including five points, since all those amplitudes are either zero, 
MHV or googly-MHV. 

It should be noted that the "Maximally" in MHV becomes a misnomer at loop 
level since the two first equations of (|2.41|) cease to hold in non-supersymmetric 
theories. 



2.3.2 Supersymmetry 

Tree-level amplitudes in Yang-Mills theory have an apparent supersymmetry 
[104,119] because, say, a gluino and any adjoint fermion have the same Feynman 
rules. Also, at the level of individual colour ordered amplitudes fundamental and 
adjoint fermions are really indistinguishable since their colour ordered Feynman 
rules are the same. The difference only shows when we compute amplitudes with 
non- adjacent fermions where the fundamental fermion amplitude must vanish. 
The actual supersymmetry only reveals itself at loop level where the number of 
particle species becomes important. 

The tree-level supersymmetry can be exploited to constrain tree amplitudes 
by "Supersymmetric Ward Identities", SWTs [95,96], which relate amplitudes of 
different external particle content but with the same momenta and the same sum 
of helicities. Following [76], SWTs are obtained by writing the amplitude as a 
string of operators <I>j cretaing helicity eigenstates by acting on the vacuum 

(0|$i«>2---«>n|0). (2.43) 

We now introduce the spinorial supercharge Qa and contract it as Q{q,9) = 
Gq°^Qa where is a Grassmann parameter and is the Dirac spinor of some 

^The word "googly" is a cricket reference often used in twistor theory. In cricket it denotes 
a ball bowled with the opposite spin of the normal, thrown by a right-handed bowler. 
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lightlike momentum. The inclusion of 9 makes Q{q,9) a bosonic operator. Since 
Qa annihilates the vacuum, we can write 



n 

= (0| [Q{q, 9),<^i'^2 ■■■<^n] |0) = Y^mi ■ ■ ■ [Q{q, 9), <^>i]--- '^>n\0), (2.44) 

=1 



where the last expression must represent a sum over amplitudes where Q has 
exchanged one particle with another. If one computes all the commutators, one 
will get 

[Q{q,9),g+{k)]=9[kq]f+, 
[Q{q,9),f+{k)]=9{kq)g+, 
[Qiq,9),f-{k)]=9[qk]g-, 

[Q{q,9),g-ik)]=9{qk)r. (2.45) 

This permits us to prove the two first formulas of (|2.41|) . The second can be 
proven by acting on {0\gif2 9t ' • •fi'n |0)> 

= {ql)A{f{, ft,---)- {q2)A{g^,gt, -..), (2.46) 

and choosing g = 1 to eliminate the first term. Notice that the amplitudes with 
two /^'s are zero because of the requirement of conservation of fermion helicity. 

We can also extend the formula for MHV gluon amplitudes to amplitudes 
with a fermion pair by acting on an amplitude with two negative helicity gluons 
at 1 and i and a positive helicity fermion at j and the rest positive helicity gluons. 
This gives the SWI 

= (<?l)^(/f ,gr,/+) + {qi)A{g^,fr,f+) " {qj)A{g^,g^,gp, (2.47) 

and choosing q = 1 reduces this to 

.4fe,/r./;) = M^te..r.^;) = - (i,)';4'"'(„i) . (2.48) 

This can be explained as moving a half unit of negative helicity from i to j, 
thereby exchanging a factor of (li) with a factor of (Ij). 

SWI's can be extended to describe other supersymmetries. For example, M = 
2 supersymmetry would contain a scalar and thus allow (I2.48P to be extended to 
include those. The M = 4 supersymmetry algebra will be of special importance 
to us since that theory enters at both tree and loop level in several places below. 
It has a gluon g^, four fermions six real scalars Sab = —Sbai a-^d the algebra 

[Qaiq,e),g+ik)]=9[kq]ft, 

[Qaiq, 0),ftik)] = 05ab{kq)g+ + 9[kq]Sab, 

[Qa{q, 9), Sbc{k)] = 95ab{kq)ft - 0Sac{kq)ft + 0[qk]eabcdfl , 

[Qa{q,0), f,^{k)] = 95ab[qk]g~ + -9{qk)eabcdScd 

[Qa{q,e),g-{k)]=9{qk)f-. (2.49) 
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Notice that we might as weh have used Sab = 2^abcdScd- Using these, we can 
deduce SWI's for theories with several different fermions and scalars, as well as 
SWI's for = 4 amplitudes to any loop order. In particular, we can deduce the 
exact form of all tree-level MHV amplitudes in = 4 SYM by a procedure like 



2.3.3 Limit Behaviour 

We will also briefly touch on the behaviour of Yang-Mills tree amplitudes in 
different limits. In a limit where a multi-particle kinematic invariant (and thus 
a propagator) goes on-shell, the residue factorizes into two on-shell amplitudes, 
summed over possible internal states. In the pure-glue case we have 

h=± 

as — > 0. This factorization is present in most theories described by Feynman 
rules and not just Yang-Mills. 

When two-particle kinematic invariants go on-shell the situation is different 
in massless Yang-Mills theory. Because of kinematics we should really be looking 
at the limit where two colour adjacent momenta become collinear. In this limit, 
the amplitude factorizes into the amplitude with one fewer external state times a 
so-called splitting amplitude which depends on the particles going collinear and 
the internal state which is going on-shell. Again, considering only gluons we have, 
when a — > zP and 6 — > (1 — z)P, 

^(...,a^^6'^^...)^ ^(...,p-^%...)Split^^(z,a'^%6'^'). (2.51) 

hi=± 

The splitting amplitudes for gluons are 

Split*i"<=(a-,6-) =0, 
Split*i'="(a+,6+) 

Split*;^^(a+,6-) 



,Jz{l-z){ah) 
y/z{l-z){ab) 



^2 

(.tree / 



Splitl-<=(a+, b-) = - (2.52) 

Y^z(l - z)[ab] 

from which the remaining can be deduced by pariaty. Similarly, one can deduce 
splitting functions involving fermions and scalars from the corresponding MHV 
amplitudes. The one-loop limits and simgularities will be considered in section 

MM 
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2.4 Gravity as a Quantum Field Theory 



It is well known that General Relativity has serious problems when implemented 
as a quantum field theory. It has a dimensionful coupling which means that it is 
non-renormalizable. At each loop order one has to introduce new counterterms to 
cancel divergences, and since nothing prevents finite parts of those counterterms, 
we automatically introduce infinitely many coupling constants into the theory. 

Because of this, the "conventional wisdom" is that a QFT version of GR is 
non-sensical. The situation, however, is not as bad as one may think. Gravity can 
be treated as an effective field theory [131], where the theory is only renormal- 
ized to the loop order needed. This is legal because we know that an all-orders 
renormalized calculation can still be Taylor expanded around zero coupling. As 
long as the coupling constant is small, we can expect tree-level computations to 
be close to the actual result, and even at loop level the counterterms contribute 
to local interactions but not to e.g. long range corrections to the gravitational 
potential [47,78,79]. One could also remark that Yang-Mills theory in more than 
four dimensions is also non-renormalizable, a fact which neither does nor should 
deter people from working with the theory. 

Another problem plaguing quantum gravity theories is that of algebra. When 
the Einstein-Hilbert action, 

Seh = ^ [ d^x^R, (2.53) 



2^2 

is written out in terms of a perturbation around flat 

Qfiu = ilf,u + Kh^i, or curved g^u = 9% + f^h^y (2.54) 

space, it contains vertices to all orders, and even the lowest order terms suffer 
from severe congestion of indices. The three-point Feynman vertex takes up 
about half a page in condensed notation, the four point requires one or two. This 
means that only amplitudes where gravitons enter in the simplest manner can be 
calculated analytically. The only computational point we will make here is that 
of polarization tensors: Gravitons are spin 2 particles which have two on-shell 
helicity states with traceless polarization tensors e^^. It so happens that one can 
use products of gluon polarization, 

4- = 4^t- (2.55) 

With the choices of polarization vectors described in section 12.21 we may even 
choose different reference momenta for and e^. 

2.4.1 The KLT Relations 

What one can do to compute tree-level amplitudes in quantum gravity is to use 
a result from string theory by Kawai, Lewellen and Tye [99], the KLT relations. 
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These are obtained by rewriting the string computation for closed string scatter- 
ing in terms of products of open string scatterings, and by taking the a — ^ 
Umit, they make a statement about a relation between graviton and gluon scat- 
tering amplitudes. The relation is valid on-shell at tree-level in any number of 
dimensions. The first three are (setting irrelevant factors to 1) 

A^(l,2,3) =^(1,2,3)^(1,2,3) 
M{1, 2, 3, 4) = si2^(l, 2, 3, 4)^(1, 2, 4, 3) 
Mil, 2, 3, 4, 5) = si2S34^(l, 2, 3, 4, 5) J(l, 4, 3, 5, 2) 

+si3S24^(l, 3, 2, 4, 5)^(1, 4, 2, 5, 3). (2.56) 

The A^'s are graviton amplitudes and the A and A are colour ordered gauge 
amplitudes. In fact, there are many forms of these relations since an (irrelevant) 
relabeling on the gravity side gives another combination of orderings on the gauge 
side. The simplest non-trivial example is the four graviton amplitude in four 
dimensions which is 

7W(1-, 2", 3+, 4+) = si2^(l-, 2", 3+, 4+)^(l-, 2-, 4+, 3+) 

(12)3 ^12)3 



(34) [43] -^2^^ ^^^^ ^^^^ ^^^^ ^^^^ ^^^^ 

l^^'K^^i' '2.57) 



(12) (23) (34) (41) (24) (13) ' 

The KLT relations also imply that graviton states decompose into a tensor 
product of two gluon states. Since the proof of the KLT relations does not 
rely on the particular string theory in question, we may imagine that there are 
more states than gravitons on the gravity side and gluons on the gauge side and 
that the gravity states are tensor products of two gauge states. The two gauge 
theories need not be the same. If we imagine the gravity side to be a heterotic 
string theory, the one gauge theory corresponds to the right moving sector and 
the other corresponds to the left moving sector [22] . 

2.4.2 The = 8/ = 4 Relation 

At this point we can take a closer look at one of the main themes of this thesis, 
namely the relation between the maximally supersymmetric gravity and gauge 
theories in four dimensions. At tree level, the KLT relations give us a corre- 
spondence between the two theories both at the level of the states of the two 
theories (something which is independent of the loop order) and between the 
actual amplitudes. 

The spectrum (and Lagrangian) of AT = 8 supergravity [72,73] can be derived 
by dimensionally reducing 11 dimensional J\f = 1 supergravity [74] and contains 
a graviton, eight Rarita-Schwinger gravitini, 28 vectors, 56 Majorana fermions 
and 70 real scalars, a total of 256 states when helicities are taken into account. 
The = 4 multiplet has a gluon, four (Majorana or Weyl) fermions and six real 
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scalars, a total of 16 states. By the KLT relation, any state in = 8 supergravity 
can be seen as a tensor product of two = 4 states such that the sum of the 
two M = 4: state helicities add up to the M = 8 state helicity. Supergravity can 
have the same construction as we saw for super symmetric gauge theories above 
where the states are related by supersymmetry charges. In that sense, four of 
the M = 8 charges map to the charges of one of the = 4 theories and the 
remaining four to the other. 

This, however, does not account for all the symmetry on the gauge side. The 
symmetry of Yang-Mills theory at tree level is the same as = 4 SYM at loop 
level, namely the superconformal symmetry. Since the spectrum, the supersym- 
metry and the amplitudes are directly related to tree-level M = 8 supergravity, 
one might guess that there is an additional symmetry in gravity theories which 
mirrors the additional conformal symmetry in Yang-Mills. We will return to the 
relations between AA = 8 and AA = 4 at loop level in chapter [71 
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Chapter 3 

Tree Level Methods 



3.1 Earlier Methods 

Calculation of tree-level scattering amplitudes is normally done by applying the 
Feynman rules. This procedure is completely safe, in the sense that it is com- 
pletely mechanical and results are trusted. The expressions generated in this way 
are known not to be the most compact, but they also avoid unphysical singular- 
ities which may cause problems for numerical evaluation. For calculations by 
hand, however, they become impractical beyond five-particle processes because 
the number of diagrams becomes too big. Taking as an example the n + 2 point 
gluon amplitude, the number of diagrams increases factorially, as shown in the 
table below [101,112]. 



n 


2 


3 


4 


5 


6 


7 


8 


# of diagrams 


4 


25 


220 


2485 


34300 


559405 


10525900 



Table 3.1: The number of Feynman diagrams contributing to the scattering pro- 
cess gg ^ n g. Extracted from [112]. 

The problem of the number of diagrams can, to a large extent, be solved 
by implementing Feynman rules (or an equivalent recursive formulation) on a 
computer, where the running time for generating an amplitude and evaluating 
it at a specific kinematical point is not prohibitively long even for a number of 
points of the order of 15. There are several programs that can do such calculations 
[98,103,110,111,123]. 

Problems may still arise if we want to generate a large number of events with 
many final states, since reevaluating the amplitude many times is time consum- 
ing. On the other hand, deriving the amplitude analytically and subsequently 
inserting numbers is also problematic since the analytic expressions from Feyn- 
man rules are extremely large. What we would ideally like to have is some means 
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of computing the amplitude in a form which is compact and well suited for in- 
sertions of mimcrical values. 

One method for this, based on Feynman rules, is Bcrcnds-Giele recursion 
[16]. It uses colour ordering and the spinor helicity formalism and in addition 
keeps as little off-shell information as possible. In this way it takes advantage 
of gauge invariance and other on-shell conditions to simplify intermediate results 
instead of waiting till the end. Berends-Giele recursion for gluon amplitudes 
works by introducing the off-shell n-point gluon current Jn which is the (gauge 
dependent) off-shell amplitude for n-gluon scattering with n—1 legs taken on-shell 
by multiplying in polarization vectors and using fc^ = 0. By looking back into the 
way Jn was calculated, we see that the off-shell leg had to be attached to either 
a three or a four vertex. Whatever was at the other legs of that vertex would 
sum up exactly to some Jm where m <n. This allows us to write down a cubic 
recursion relation for J^, which can be terminated by setting the momentum of 
the off-shell leg on-shell and multiplying by the polarization. We can either try 
to solve the recursion relation exactly or implement it on a computer. A good 
description of Berends-Giele recursion can be found in [76] . 

Even though this sort of recursion can speed up the calculation of ampli- 
tudes, it still suffers from very long expressions, although less severely than the 
traditional application of Feynman rules. 2004 saw the invention of two methods 
which, in their separate ways, have contributed to both more compact results 
and theoretical insight into the structure of gauge theory amplitudes. We now 
turn to them in historical order. 

3.2 MHV Rules 
3.2.1 Bcisic Construction 

As described in the introduction, the analysis of Yang-Mills amplitudes in twistor 
space led to an intuitive picture of MHV vertices connected by propagators. 
The MHV vertex, however, was not a well defined concept since it is only well 
defined on-shell where the momenta are all light-like. To make the MHV picture 
concrete, it was necessary to construct an off-shell continuation of the MHV 
amplitude. This was achieved by Cachazo, Svrcek, and Witten [69] in March of 
2004. The rules presented below were not derived as such, but rather proposed 
and subsequently motivated. For multi-gluon amplitudes, the rules state: 

1. Draw all possible graphs ("MHV diagrams") 

• where external lines represent external gluons and have their helicity 
marked, 

• where internal lines have opposite helicities (it) at each end, 

• where each vertex has exactly two edges of negative helicity — and at 
least one of positive helicity — attached. 
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• and which is planar and respects the colour ordering of the external 
gluons. 

2. For each graph, assign 

• to each external line the outgoing momentum of the respective gluon, 

• off-shell momenta Pi to the internal lines by assuming conservation of 

momentum at each vertex, 

• holomorphic spinors to internal lines as 

i^^'') = im (3.1) 

where the same (almost arbitrary) anti-holomorphic spinor |j7] must 
be used for all internal lines of all graphs. 

3. For each graph, multiply, 

• for each internal line of momentum Pj, 

1 

^' 

• for each vertex, the MHV formula using 
spinors of the internal lines. 

4. Add the contributions of all graphs. 

The combinatorics of these rules are similar to that of Feynman rules. \mt 
with a somewhat gentler sprawl of diagrams. It is simple to show that an MHV 
diagram with p negative helicity external gluons will have p — 1 MHV vertices 
and p — 2 propagators; significantly fewer than for Feynman rules. For fixed p, 
the number of diagrams grows polynomially in the number of external gluons, 
thus presenting a significant simplification over Feynman rules, but if we need all 
helicity configurations (and thus all p up to half the number of gluons) e.g. for 
computing an unpolarized cross section, the asymptotic behaviour of the number 
of MHV diagrams is again factorial, thus not providing a significant advantage 
over Feynman rules for many external gluons. 

For few gluons, however, the simplification is remarkable. By taking advan- 
tage of parity, a six gluon amplitude can be computed with the use of only 32 
diagrams (depending on how you count) three of which are MHV and thus trivial. 

To gain some familiarity with the rules, let us consider the simplest non- 
trivial case of A(l~, 2~, 3~, 4+, 5"*"). It is of course slightly trivial because it is 
a googly-MHV amplitude, but it will be a welcome first check to verify that the 
rules obey parity. Having three negative helicity gluons, this amplitude is called 
"Next-to-MHV" , or NMHV, and the diagrams must have two vertices and one 
propagator. There are four permitted diagrams. 



(3.2) 

\P^) for the holomorphic 
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Following the rules above, we get the contribution of the first diagram, 

((-P5^)5)(5f)P|i(34)(4P5^)(P5^2) 
{IP51V? 1 (23)3 



(5P5i??](51) [15] (5f) (34)(4P5ir?](2P5ir/] 

[5r]]^{23f 

[lr/][51](34)(4(5 + 1)77] (2(5 + l)r/]- 



(3.3) 
(3.4) 
(3.5) 



At this point we can make a choice of \t]] if we like. The choice \r]] = |5] will 
immediately set the diagram to zero, where as the choices \r]] = |1], |(5 + 1)4), 
1(5 + 1)2) will render it undefined. The latter three correspond to unphysical 
poles which cancel among the diagrams. If we continue with |r/] = |5], the three 
other diagrams contribute 

(34)^[45]3 ^3^^^ 



[12][51](45)[52](5Pi2 5]' 

(14)^[45]3 

[23] [51] (45) [52] (15) [53] 

and, 

_ (12)2 [45]3 

[34][51](15)[53](5P34 5]' 

which adds up to the expected result, 

[12] [23] [34] [51]- ^ • 



(3.7) 



(3.8) 
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Though the result of this calculation was already known to be fairly simple, 
it is remarkable that we can compute a five point amplitude using only three 
contributing diagrams of such limited complexity. 

A comment on the choice of the off-shell extension (j3.ip is in order here. 
Because internal lines must have opposite helicities at each end, every diagram is 
invariant under scaling of internal spinors. In other words, the requirement (13. ip 
should really read 

\Pt) ~ |P,r?]. (3.10) 
As noted by Kosower [102] such a relation can be achieved by defining 

p\^' = Pl' + at (3.11) 
and choosing a such that is lightlike. This results in 

p^, ^p,_ (3.12) 

\Pt)[Pl\J-^^^ {3.13) 

In the Yang-Mills MHV rules, only the holomorphic spinor is used, so in a sense, 
the above rewriting adds information which is not used in practise. 



3.2.2 Gravity: A Puzzling Failure 

The success of MHV rules for Yang-Mills theory immediately raised the question 
of whether MHV rules existed for gravity. Even though such rules were not di- 
rectly implied by the KLT relationship, the possibility was not excluded by them, 
nor by any other known principle. The main obstacle, however, was that gravi- 
ton MHV amplitudes were not holomorphic (c/. p. 560 ) and thus the extension 
\P^) ~ l^^i^] was insufficient. As an example, take the graviton equivalent of the 
diagram treated in detail above. MHV rules would give that this diagram was 

M(5+,l-,-P5\-)-^^M(2-,3-,4+,P5^+), (3.14) 

where the graviton amplitudes can be obtained from the KLT relations. For that 
we would need ^2^3 o^") equivalently, (P4 + P^i)'^, but these are not obviously the 
same because we are now "inside the MHV vertex" where we must use rather 
than P51. For this to be consistent, we must have 



(4P5^4] = (4P5i4] + P|i 


(3.15) 


(4P5ir?](r?P5i4] _ 2 


(3.16) 




(3.17) 


{VP451V] = 


(3.18) 


(r?(2 + 3)??] = 0. 


(3.19) 
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The last requirement states that rj^ must be a hghthke momentum orthogonal to 
the time-like momentum P23. By going to the rest frame of P23 it can be seen 
clearly that no such (non-zero) ry^ exists. Thus, any gravitational attempt at 
MHV rules along these lines seems impossible. 

Although it seems that the problem at hand is rooted in the non-holomorphic 
character of graviton MHV amplitudes, it is in fact much deeper. The amplitude 
M{2~ ,3^ ,4'^ , P^f) is ill-defined because we have taken P51 away from its real 
value without compensating in other places; the amplitude isn't just off-shell, it 
disobeys conservation of momentum. Had this been an amplitude proper, the 
implicit delta function in the momenta would have set it to zero before we even 
started. And this argument relates as much to Yang-Mills as to gravity. 

In this light, the fact that Yang-Mills MHV rules do work becomes even 
more puzzling. It gives the impression of being an inherently ill-defined proce- 
dure which magically becomes meaningful for theories with holomorphic MHV 
amplitudes. The missing ingredient to understand why the situation is signifi- 
cantly brighter than this, will appear in section [3^ and the puzzle will be resolved 
in the next chapter. 



3.2.3 Fermions, Higgses, and Vector Bosons 

Meanwhile, we can turn to some of the additional successes of MHV rules. As 
described above, the success of the MHV rules seemed to rely somewhat on the 
holomorphic nature of MHV amplitudes, so similar theories with holomorphic 
MHV amplitudes ought to have associated MHV rules. 

The primary example of this is Yang-Mills theory with coloured (adjoint 
or fundamental) massless scalars and fermions [90,91,134,135]. As described in 
section [2.3.2t replacing gluons with fermions or scalars in an MHV amplitude only 
changes the numerator, and there is nothing preventing us from extending the 
MHV rules with fermionic and scalar internal lines. Particles in the fundamental 
representation are handled by requiring them to be adjacent at all times. If, for 
instance, we wanted to compute the NMHV amplitude A^{q^ , 2^, 3~, 4+, g^), the 
diagrams would be similar to those of our previous example, just with some of 
the gluons being replaced by fermions denoted by dashed lines. 
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4+ 




> 



+ 




4+ 



> 



+ 




3" 



The contribution of the second diagram would be 



(12)3((-P|2)2) 1 {3Pl,f{35) 



(3.20) 



(12)2 1 (3^^2)2(35) 



(3.21) 



(lPt2)Pf2(34)(45)(5P|2)' 



which differs from the corresponding gluon diagram by the factor 



(2P,^,)(35) _ [lr/](35) 
(21)(3P|2) {3Pi2V]' 



(3.22) 



If we now try to compute ^45(51, 2 ,3 ,4+, 55) we will again have four diagrams 
similar to those above and for the one treated we would obtain the factor 



relative to the gluon result. 

Another case where MHV rules apply is one to which we will turn in chapter [SI 
namely that of a massive, uncoloured particle H (such as a Higgs) which couples 
to a gauge field through the term CHF^^F^^ in the Lagrangian. 

The last case is that of the coupling of gauge fields to electroweak vector 
boson currents. The formalism is somewhat involved, and rather than trying to 
explain it here, we refer the reader to the original article on the matter [44]. 

3.3 On-Shell Recursion Relations 

As was discussed in the introduction, a decisive turn came with the realization 
that spinors could be taken away from 3+1 dimensions, in which Aa and are re- 
lated by complex conjugation, to complexified Minkowski space where they are 
complex and independent. This allows for the use of standard methods of complex 
analysis, in particular Cauchy's Theorem. Although on-shell recursion relations 
were first discovered as the combined consequence of = 4 IR consistency con- 
ditions and the quadruple cut method [60], we will take the complex analysis 
viewpoint introduced by Britto, Cachazo, Feng, and Witten (BCFW) [61]. We 
will, however, also touch on the connection with IR consistency conditions in 
later chapters. 




(3.23) 
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3.3.1 Using Cauchy's Theorem on Yang— Mills Amplitudes 



BCFW on-shell recursion works by deforming a tree-level amplitude ^ by a com- 
plex parameter z to A(^z) by choosing two external particles i and j and deforming 
their spinors 

l^] = K] + ^|j], \i) = \i)-z\i). (3.24) 

This deformations obeys conservation of momentum since 

Vi^Vj = \i)[A + = + z\i)[j\ + - z\i)[j\ =pi+pj, (3.25) 

so A{z) is a well defined amplitude in C^. 

For now, we assume that A{z) as z ^ oo. This allows us to use Cauchy's 
Theorem to write 



ReSp^(2;) 



(3.26) 

poles p ^ 

The poles of tree amplitudes are well known; they correspond to the Feynman 
propagators of internal momenta going on-shell, and the "residue" is given by 



lim 

P? ^0 



Ak, . . . , m, -P^X^A{P^^^, m + l,...,k- 1).(3.27) 



h=±l 



If i G {k, . . . , m} and j ^ {A;, . . . , m} (or vice versa) there will be a corresponding 
pole in z because 

Pk,rn{zf = {Pk,rn + z\i)[j\f = + z{iPk,mj], 

thus 



p2 

k,m 



{iPk,m3] ' 

The residue of this pole is 
zA 



(3.28) 
(3.29) 



lim 
1 1" 



E 



Aik, ...,m, -P'kX(zk,m))AiP^^Jzk,m),m + l,...,k-l) 



h=±l 

Combining this with (I3.26|) leads to 



{'iPk,mj] 



(3.30) 



A 



E 



Aik, . . . , m, -P^^{zk,m))A{P^^{zk,m).m + 1, . . . , /c - 1) 



e {i + 1, . . . , i} 

m g {i, . . . , j - 1} 
k ^ m,m + 2 
h = ±l 



p2 

k,m 



(3.31) 
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In human language, the conditions under the sum state that we must sum over 
all internal momenta affected by the deformation. We can draw this as a diagram 
consisting of two deformed amplitudes and an internal propagator, together with 
the sum, 



E 

all such diagrams 



, ± 



Pi 



k,m 



Pk,m is on-shell and has the form 

Very often this recursion relation is used with i and j adjacent, e.g. = (n, 1) 
in an n-point amplitude. This simplifies the above relation to 

, ^ A{m+l,...,n, P{^^{z^))A{-Pl^{zm), 1, . . . , m) 

•A-n - 2^ 2 {6.66) 

m=2 h=±l 1."^ 

where 

^l,m 3 \Pl,mn]{lPl^m\ \Pn,mn]{lP2,m\ 



Before we can go and apply this decomposition on poles to the calculation 
of real gluon amplitudes, there are two subtleties to consider, having to do with 
three-point amplitudes and the behaviour of ^(z) as z ^ oo. 

3.3.2 Three-Point Amplitudes (and an Example) 

In real Minkowski space, three-point gauge amplitudes are zero. This basically 
stems from the fact that in order for the three momenta to be light-like, they 
must also be orthogonal, 

= pi = (-pi-P2)' = (12)[21], (12) = [12]* = 0, etc. (3.35) 

=^-pi~P2~P3- (3.36) 

When we move away from real Minkowski space, however, the restriction that 
Aa = Aq is lifted and the on-shell constraints have two solutions, 

|1) ~ |2) ~ |3) or |1] ~ |2] ~ |3]. (3.37) 
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These each ahow one of the two hehcity configurations to be non-zero because 
the expressions 

[12]3 (12)3 

and . , (3.38) 



[23] [31] (23) (31) 

respectively, are weh-defined and contribute to ()3.3ip . 

We can illustrate the method and the considerations about three-point ampli- 
tudes with the calculation of the googly-MHV five-point amplitude ^(1~ , 2+ , 3~ , 4''" , 
We define the deformation by taking 

|3] = |3] + z|4], |4) = |4)-z|3). (3.39) 

The internal momenta affected by this deformation are Pi^a and P2,3 each corre- 
sponding to a pole in z at 

^',3 _(45) _ Ph_ [23] 

'''' - - (37^ - (35) ' '^'^ " ~W^] ~ ~WY ^ ^ ^ 

Summing also over internal helicities gives four diagrams 





with the contributions 

^(l-,2+,3-(zi,3),-Pit3(^l,3))^(7^i:3(^l,3),4+K3),5-) 



p2 

^1,3 



■4(l',2+,3-(zi,3),-Pi73(zi,3))-4(P+3(^i,3),4+(zi,3),5-) 

M,3 

■4(2+, 3-(z2,3), -P2t3(^2,3))^(\3K3),4+(z2,3), 5', 1") 

^2,3 

■4(2+, 3-(z2,3), -P2:3(^2.3))^(^2"!"3K3),4+(Z2.3), 5', 1") 

p2 
^2,3 



(3.41) 
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Two of these contributions can be discarded immediately. Number 2 and 3 con- 
tain an amplitude of the form A{ h) which is zero. In contribution number 

1, zi^3 is found by requiring that 

A% = (45)[54] =0, or |4) ~ |5), (3.42) 

so the three-point MHV amplitude in that contribution is zero, even in complex 
momenta. The corresponding condition in the last contribution is |3] ~ |2] which 
kills the googly-MHV three-point but not the MHV three-point we need. Thus 
it gives the only contribution, for which we need the solutions 

, V [231, , IP342I - |i^2 34l(3P23l l^'2 44l[2| 

This allows us to compute 

(3(-P2,3))^ 1 (51)3 



((-P2,3)2)(23) Pis (1P2,3)(^2,34)(45) 

(51)=^(3P2,3)=^ 
(23)2[32](45)(2P2,3)(1^^2,3)(4P2,3> 

(51)3(3P2,34]3[24]^ 

(23)2 [32] [2P3,45) (2P2,34] (1^2,34] [2P3,4P2,34] 

(51)=^(32)3[24]5 

(23)2 [32] [21] (15) (23) [34] (15) [54] [24] (51) [15] 
[24] 4 



(3.44) 



[12] [23] [34] [45] [51]' 

which is, of course, the known expression for a five-point googly-MHV amplitude. 



3.3.3 Asymptotic Behaviour of the Deformed Amphtude 

Until now we postponed the issue of whether lim^j^oo -^{z) = 0. For certain 
choices of the spinors to deform we can immediately argue that this must be 
the case by looking at the calculation from Feynman diagrams and tracking the 
dependence on z. Every propagator in a Feynman diagram affected by the defor- 
mation must asymptotically contribute z"^, and every three vertex affected by the 
deformation must contribute z^. The "worst case scenario" for the asymptotic 
z dependence will be when the deformation only affects three- vertices. Before 
multiplying by the external polarizations, such a diagram will go as z^. 

If we make the deformation ()3.24p the dependences of and on z will 
depend on their helicities. For i the reaction to the deformation is 

e^^ = ^^z\ e^~ = ^^z-\ (3.45) 
V2{iq) ' V2M ^ ' 
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and for j the asymptotic behaviour is the opposite. This leads us to conclude 
that the worst behaviour of A{z) as z ^ oo for helicities {hi, hj) is 



-1 



(+,+),(-,-) ~z 



1 



(+,-) 



~ z' 



,3 



(3.46) 



This tells us that using BCFW on-shell recursion relations is always justified 
when we deform the antiholomorphic spinor of a negative helicity particle and 
the holomorphic spinor of a positive helicity particle. 

In practice, however, the behaviour may be more benign. If we consider the 
MHV amplitude A{1~,2~, . . . ,n"'"), the choice (i,j) = (1~,2~) goes as z~^ and 
so does e.g. the choice (3"*", 4+) while (2~,4"'") and (3^,5+) go as This shows 
us that we can often expect an actual behaviour which is a few powers of z better 
than the worst case. This is another expression of the fact that Feynman rules 
are not taking advantage of all the available symmetries of the theory. 

3.3.4 Some Consequences of BCFW on-shell recursion 

Two points are worth making at this point. First, since we have rigorously 
proven that we can always perform on-shell recursion with the choice (— ,+) it 
follows that any tree amplitude can, by successive recursion, be reduced to three- 
amplitudes. In this perspective, the four gluon vertex can be completely neglected 
for practical calculations, and its presence in the Lagrangian is only required for 
gauge invar iance. 

Second, the existence of on-shell recursion shows that other formalisms that 
claim to reproduce Yang-Mills amplitudes are correct if they have the right sin- 
gularity structure (13.270 and Lorentz covariance. This was proven together with 
the proposal of the MHV rules [69], thus, BCFW recursion provides an indirect 
proof that the MHV rules are correct. 

Both of these points are elaborated on in the original article [61]. 

3.3.5 Extension to Other Theories and Deformations 

The deformation (13.240 chosen by BCFW is by no means unique. For instance, 
one may choose three gluons i,j,k and deform them to [39] 



which also obeys momentum conservation for any a. In general one will have to 
redo the derivations of section 13.3. H but most features of the relations persist. 
Using extended deformations may be helpful for improved z ^ oo behaviour and 
for constructing relations to the MHV rules as done in the next chapter. 

The methods of deriving on-shell recursion relations are clearly very generic, 
relying mostly on the fact that the theory at hand is an ordinary quantum field 
theory. Gauge invariance or colour ordering may be practical tools for obtaining 
simpler expressions, but are not required as such for doing on-shell recursion 
in this way. One extension is to allow other states than gluons to participate. 





(3.47) 
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such as charged scalars and fermions (massless, massive, D-dimensional) [7, 9, 
107,108,125], another is to take a completely different theory, such as QED [117], 
gravity [12,68], or more exotic theories [15]. The challenge in most of these cases 
is to prove the appropriate z — > oo behaviour, which is sometimes done and 
sometimes left to faith. Especially in gravity, where the "worst case scenario" 
analysis gives very bad asymptotic behaviour, the on-shell recursion relations 
were used for quite some time on this faith until a proof appeared [14]. 

The applications to loop-level calculations will be reviewed in section 15.51 
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Chapter 4 

MHV Rules from Recursion 



In the present chapter we will connect the two subjects of the previous chapter and 
show how MHV rule constructions follow from on-shell recursion. The original 
articles are [126] for the Yang-Mills case and [49] for the gravity case. 



4.1 MHV Rules in Light of Recursion Relations 



The MHV rules as described in section 13.21 distance themselves from recursion 
relations in one obvious way: recursion relations have one unshifted propagator, 
while the MHV rules have several. This immediately tells us that a simple deriva- 
tion of the MHV rules from recursion is not obtainable, but also that the NMHV 
case can serve as an initial test case. With this in mind, we can note four features 
of MHV rules that will guide us: 

1. An internal momentum is shifted on-shell in a manner close to that of 
recursion relations. Setting 

P' = P-:r^V, (4.1) 

is practically the same as shifting P as 

P = P- zr], (4.2) 

setting P^ = 0, and solving for z. 

2. The asymmetry between holomorphic and anti-holomorphic spinors is deeply 
rooted. The attempt to construct MHV rules for gravity by choosing |r/) = 
{\rf\)* failed because it produced non-Lorentz invariant results, thereby pro- 
viding a hint that r/ is a complex momentum whose holomorphic and anti- 
holomorphic spinors are unrelated. Admittedly, there is a bit of hindsight 
in this argument, since the existence of MHV rules for gravity was not 
guaranteed. 
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3. If 7] is a complex momentum, we might as well continue all momenta. If 
we do that, all the anti-holomorphic spinors of the external momenta are 
unspecified by the MHV rules, just as \r]). In fact, \r]) may even differ 
between different internal momenta. 

4. The MHV rules contain MHV three-point vertices and no googly three- 
point vertices. The easiest way of ensuring this would be to shift only 
anti-holomorphic spinors. This fits well the fact that these do not show up 
in the MHV rules anyway. 

This leads us to conclude that if MHV rules can be proven from recursion re- 
lations, the proof must involve a non-minimal shift of anti-holomorphic spinors, 
and all these shifts must be proportional to the (otherwise unspecified) \rj]. 

4.2 The NMHV Case 

We first apply the above considerations to the calculation of the amplitude 
Anini^ , . . . , 1712 , • • • ' '^s ' • • •) where the . . . have helicity -|-. If we choose to shift 
only the three negative helicity gluons, the shift is unique up to a constant. 



Imi] = \mi] + z\ri]{m2ms) , (4.3) 

1^-2] = \m2] + z\r]]{m3mi) , (4.4) 

1^-3] = Ims] + z\r]]{mim2) , (4.5) 
where conservation of momentum follows from the Schouten identity 

|mi)(m2m3) -|- |m2)(?n,3mi) -|- \m3){mim2) = 0. (4.6) 



We will postpone the proof of the legality of this shift to section 14.51 With the 
above shift, the z dependent internal momenta are those that have at least one 
of ?7ii,2,3 at each end, e.g. 




The opposite helicity assignment for the inner line does not contribute, so we 
have split the amplitude into two (shifted) MHV amplitudes. The solution for z 
is found by setting 

= {P — z\mi) {m2m^)[ri\)'^ = P^ — z{m2m^){miPr]] = (4.7) 
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(4.8) 



It should now be clear that we have reproduced the NMHV version of the MHV 
rules for Yang-Mills: We must sum over all diagrams that split the amplitude 
into two MHV vertices. For the external momenta we must use the shifted spinors 
(but because the MHV amplitude only involves holomorphic spinors, it makes no 
difference) and for the internal we must use |P) ~ \Prj\. 

We have, however, obtained more than just the Yang-Mills MHV rules. We 
have obtained the full expression for the shifted momenta, allowing us to con- 
struct MHV rules for theories with other MHV amplitudes than Yang-Mills; most 
notably gravity. We will return to the MHV rules for gravity in section 14.41 after 
treating the general case with more than three negative helicity particles. 



4.3 The NNMHV and General Cases 

From the above discussion, it should be clear that setting up a recursion by 
shifting the anti-holomorphic spinors of all negative helicity gluons by the same 
reference spinor \rj\ allows a description where a) the left and right amplitudes 
both have fewer negative helicity gluons than the full amplitude, b) there are no 
three-point googly amplitudes, and c) the spinor to be used for the internal line 
is the same as in the MHV rules. The remaining problem is that there are several 
propagators in each MHV diagram while the recursion relations can only provide 
one. 

The solution will be to perform similar shifts one after the other and to reach 
a reperesentation which can subsequently be proven to be equivalent to the MHV 
rules. Starting with an NNMHV amplitude, we can choose a set and shift 

\fhi\ = |mj] + zcf^rj]. (4.9) 
For conservation of momentum, the c's must be chosen such that 

4 

^cf\mi\=Q. (4.10) 
1=1 

If we label the possible internal lines by a, we can write the result of the recursion 
as 

^„(,„- J . •^'"'"•^"'^f (4.11) 

a » 

where \Pa) = \Paff\-i ^^^d there is an implicit summation over helicities where 
relevant. The two sub amplitudes are MHV and NMHV, respectively. Such 
contributions can be depicted as 
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To uncover a second propagator, we can set up a similar recursion for each term 
in the above sum. However, we have the constraint on those recursions that 
we should not uncover the same propag ator 

again, or we would not have 

(2) 

progressed. Thus, for each a we choose a set ■ and set up the recursion as 

[m,|+zcg[r?| (4.12) 



nii 



such that P2 = p2 This gives one additional constraint on the c'^j apart from 
conservation of momentum. Performing the recursion for all terms in the sum 
()4.1ip . we arrive at 



P'aPi 



(4.13) 



where the terms in the sum can be drawn as (placement of m's and it's can differ) 




Since our recursion cannot split MHV amplitudes, it is clear that all three ampli- 
tudes in the numerator must be MHV and thus that the helicities of the internal 
lines are fixed by the external helicities. Moreover, if we write the twice shifted 
|m]'s as 



\fhi] = \mi] + ri\ri\ 



(4.14) 



and remind ourselves that the coefficients were constructed to obey momentum 

^2 ^2 ^ 
conservation and = = we can fix the \ fhi\ uniquely as a function of a 

and h. This solution is independent of which of these was uncovered first, so we 

can collect terms in (j4.13p two by two and write it as 



{a,fe} 



-p. 



(4.15) 
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1 1 

XI + 



p2 p2 p2 p2 

.(1) 



(4.16) 



where Pb{a) means Ph shifted using c\ as to put Pa on-sheh. The last factor, 
however, is exactly what we would get if we took 



p2p2 



(4.17) 



and did recursion on it using ci \ In terms of diagrams, what we are arguing is 
that 




+ 




m4 



m.4 




m4 



Thus, we can write the amplitude as 

(m 



A.(.nr 4) = E -^'^"■^'>-^<^'--fP>')-^<^"--P>') . (4.18) 

{a,b} ^"^b 



This is exactly the MHV rules for NNMHV amplitudes. 

It should now be clear how to proceed with the general case. To uncover the 
first propagator in an N^MHV amplitude, a, use the shift 

\mi] = \mi\ + zc\^'>\ri]. (4.19) 

For each a, uncover the next propagator b by using the shift 

Im] = \mi] + zci^^Jrj] (4.20) 

(2) ^2 

where c-^ is chosen such that P„ = 0. The third propagator c can now be 

uncovered by using the shift 

\m] = \m] + cf'l,\v] (4.21) 



45 



^2 ^2 

where c^^j^ is chosen such that = Pj, = 0. We can evidently continue in 
this way until we have uncovered p propagators, at which point we have shifted 
the \mi] to \mi] + ri\ri] and required that the shifted momenta of the uncovered 
propagators are on-shell. The latter, together with conservation of momentum, 
fixes Tj independently of the order in which the propagators were uncovered. 
Writing the amplitude formally as 

/I MHV J MHV yiMHV 

An= J2 ^^^—^ ' (4-22) 

ordered P^P^P^ ■ ■ ■ 

sets of p 
propagators 

this allows the formal rewriting 

An= J2 ^MHV^MHV...^MHV J2 . (4.23) 

unordered orderings P^P^P^ - ■ ■ 

sets of p 
propagators 

The last sum is exactly what one would get by considering 

^ (4.24) 



PiP^P^ ■ ■ ■ 



and applying recursion, first using Cj^\ then, depending on which propagator 

(2) 

was uncovered, applying recursion using q ^ and then (depending on a and b) 

(3) ' 

applying recursion using c- ^ ^, etc. In other words, 

/I MHV /I MHV /I MHV 
An= E p1p2p2; ^ (4-25) 

unordered a b c 

sets of p 
propagators 

We can now state what we have derived: To compute an amplitude, 

1. Draw all permitted MHV diagrams, 

2. For each diagram, shift the negative helicity particles according to 



|m,-l = Irrii 



+ ri\ri] (4.26) 



and solve for under the condition that all shifted internal momenta P„ 
are on-shell, and under the condition of conservation of momentum, 

J2ri{mi\=0. (4.27) 
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3. For each MHV vertex, multiply by the amplitude as a function of the shifted 
momenta. 

4. For each internal line with unshifted momentum Pa, multiply by l/Pa- 

For pure Yang-Mills, step 2 is irrelevant because it only provides information on 
the anti-holomorphic spinors. For internal momenta, \P^) ~ \Pt]] is independent 
of the actual values of the rj. In the gravity application of these rules, to which 
we now turn, this step becomes crucial. 



4.4 MHV Rules for Gravity 

Postponing again the question of the required large z behaviour, we can apply 
the above version of the MHV rules to gravity amplitudes. We will do so through 
the example of the amplitude M6(l~, 2", 3", 4", 5+, 6+), which is N^MHV, but 
also googly MHV, a fact which will not concern us here. This amplitude has a 
total of 90 MHV diagrams that are permutations of six basic ones, 



+ 




To compute the contribution of a representative of the first basic diagram. 
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1- 



6+ 



we first write down the conditions for the on-shell'ness of the shifted internal 
momenta, 

{P^f = {Pa + ri|l)[r?| + r2\2Mf = + ri(lP„r?] + r2(2P,r?] 

= (12)([21] + ri[27?] + r2[r]l]) = 0, (4.28) 
= (P, + r4|4)[77|)2 = p2 + ^^^4(5 + 6)r/] = 0. (4.29) 

Two additional conditions can be obtained by hitting the conservation of mo- 
mentum equation with two spinors (here we choose 3 and 4): 

n(13) + r2(23) + r4(43) = (4.30) 
ri(14) + r2(24)+r3(34) =0. (4.31) 

The solution is 

(4(5 + 6)1] , \{5 + 6)A)m 

(4(5 + 6)r?]' [??(5 + 6)4) ' ^^'""^^ 

^ (4(5 + 6)2] , ^ l(5 + 6)4)[r?2] 

' (4(5 + 6)r?]' [r?(5 + 6)4) ' ^^'"'■'^ 

(4(5 + 6)3] , |(5 + 6)4)[r?3] 

^ (4(5 + 6)r/]' '""J [ry(5 + 6)4) ' 

r. = _ (i±i±^ |4h = |(5 + 6)(4 + 5 + 6)r?] 
" (4(5 + 6)r?] ' J (4(5 + 6)r/] ' ^ ""^^ 

such that 

pb _ |(l + 2)r?](4(5 + 6)| 

(4(5 + 6)??] ' ^ ""^^ 

, _ |(4 + 5 + 6)7?](4(5 + 6)| 

' - (4(5 + 6)r?] ■ ^^-^^^ 

The contribution of this diagram is thus 

^^^MHV(pb-^3b-^pb+) 1 ^MHV(_p^b-^4b-^5+^g+) (4.38) 

-'ft 



(12)6 1 ^p^b3^6 1 [56](P,^4)6 
(2P,^)2(P^1)2 P2 (3pb)2(p^bpb)2 p2 (56)(45)(6P,^)(46)(5P,^) 



(4.39) 
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(12)^ 1 [r?(l + 2 + 3)3)6 

[lr?]2[r/2]2 (12) [21] (3(4 + 5 + 6)r?]2[r/(4 + 5 + 6)(1 + 2)r/]2 

1 [56] [?7(5 + 6)4)6 

^ (4 + 5 + 6)2 (45) (56) (46) (6(4 + 6)77] (5(4 + 6)??] 

(12)[56](3Pbr/]2(4P57?]6 

[lr?]2 [2r?]2 [3r/]2 [21] (45) (56) (46) {QPi^r,] {5Pb7]]P^ " 



(4.40) 
(4.41) 



Though the procedure to obtain one out of six contributions to the amphtude is in 
itself quite comphcated, it is nothing in comparison to the Lagrangian approach. 
But again, the on-shell recursion relations for gravity will produce results faster. 

Rather than computational speed, the interesting thing here is that it is pos- 
sible to formulate MHV rules for a gravity theory. This highlights one of the mys- 
terious connections between Yang-Mills theory and gravity which is not properly 
understood. In this case, the MHV rules for YM drop out of a description in 
twistor space, and such a description does not seem likely for gravity, although 
there have been attempts which reproduce the correct spectrum [1] but not the 
correct dynamical theory [116]. Thus, an explanation outside of on-shell recursion 
relations is still an outstanding problem. 

4.5 Large z Properties of the MHV Rule Shift 

The question of the z ^ 00 behaviour of deformed amplitudes should now be 
addressed. To see what behaviour we need, consider the very last step of the 

successive recursions where, for an N^MHV amplitude, we have uncovered n — 
1 propagators. To uncover the last propagator, the last shift must have the 
amplitude go as z~^. When wc get to this point, the antiholomorphic spinors 
have already been shifted many times, 

Im] Im] + {si + tiz)\r]] (4.42) 

where Si are the consequences of all the previous shifts and the U impose con- 
servation of momentum and the vanishing of the z dependence of all propagator 
terms we have already uncovered. Had we not imposed this, every propagator al- 
ready uncovered would have contributed to the large z limit. In other words, 
had the ti been unconstrained (apart from momentum conservation) the large z 
behaviour of the final shift would be z~'^. The dependence on the Si is sublead- 
ing in z and mainly a consequence of the c\ ^ 's chosen along the way, so their 
presence should not disturb this conclusion for generic external momenta. This 
leads us to conclude that the condition for the MHV rules to hold as described 
here is that under the shift 

Im] \m\i + iiz\i]] (4.43) 

where ti impose momentum conservation but are otherwise free, an N^MHV 
amplitude must go as z"" as z 00. 
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For gluon amplitudes in Yang-Mills this can be proven directly by looking 
at the worst possible Feynman diagram. In such a diagram, the n + 2 external 
negative helicity gluons must be connected by some number v of three-vertices 
and V — 1 propagators. Since three- vertices go as and propagators go as z~^, 
these internal parts of the diagram go as . The dependence of the polarizations 
is each, because we are shifting the antiholomorphic spinor of a negative 
helicity gluon, so the polarizations together go as z~^^^'^\ In total the whole 
amplitude goes as 2:~("'+^) which is one power better than required. This forms 
the last ingredient of the on-shell recursive proof of the MHV rules for gluons. 

For gravitons it is a different story. If we do the same analysis as above, 
we find that the v vertices contribute z^", the propagators contribute z~^^^ and 
the external polarizations contribute z~'^^'^~^'^\ This gives in total ^i^"^""^. The 
maximal number of vertices that may have z dependence in an m-point amplitude 
is m—2, so the worst possible Feynman diagram has the behaviour ^"^"^""^ where 
we require z"". An analysis using the above result for Yang-Mills together with 
the KLT relations gives the same estimate. 

Luckily, these types of arguments are known to be way off for gravity. If 
we do the same analysis for the BCFW shift (|3.24|) of an m-point amplitude we 
would obtain z"^~^, where the real result is known to be [5,14]. Cancellations 
of this order of magnitude in the shifts used here should in general ensure that 
the MHV rules for gravity are justified. In the original article on gravity MHV 
rules [49], the NMHV six and seven point amplitudes were found to go as z~^ 
which is two powers better than needed, and four and five powers, respectively, 
better than the worst expectation. If we allow ourselves some speculation, this 
hints that the actual z — > oo behaviour is z"""^, implying that m — n + l powers 
are removed by obscure symmetries in the gravity theory. 
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Chapter 5 

Loop-Level Methods 



Until now we have mostly concerned ourselves with tree-level amplitudes. For 
most practical applications though, amplitudes at one or more loops are required 
for interesting results. This chapter will review four general tools out of strict 
historical order. For loop-level calculations it is often necessary (or just useful) 
to use combinations of them. 

5.1 Structure of One-Loop Amplitudes 

First, however, we will start out by reviewing what is known in general about 
one-loop amplitudes, mostly in Yang-Mills theory, but also more generally. 

5.1.1 Integral Reduction 

When doing one-loop calculations with Feynman rules, one will encounter inte- 
grals with arbitrary numbers of propagators and tensors depending on the loop 
momentum of arbitrary rank. In Yang-Mills theory, an n-point calculation will 
give up to n loop propagators and up to rank n tensors, while e.g. a gravity 
calculation will give Tip to n loop propagators and rank 2n tensors in the loop 
momentum. Such integrals are impossible to do directly, so one will normally do 
some sort of integral reduction. 

The most important of these is Passarino-Veltman reduction [121], which 
reduces the powers of loop momentum appearing in the numerator by using the 
knowledge that the result must be writable as a Lorcntz tensor, together with 
the knowledge of which Lorentz tensors it can depend on. As a simple example, 
take the integral 



Since this integral can only depend on g^y and P^, it must be writable as 




(5.1) 



A{P)g^''' + B{P)P^Pr 



(5.2) 
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By contracting with g^^, and on both sides of this equation and using similar 
results from simpler integrals, one can find expressions for A and B in terms of 
integrals fewer powers of loop momentum in the numerator and the same number 
of, or fewer, loop propagators. Iterating this procedure can remove all dependence 
on the loop momentum in the numerator and leave us with scalar integrals only. 

Depending on the dimension of spacetime, the scalar integrals may be reduced 
further. Take for instance the massless hexagon integral in four dimensions (with 
implicit ze's) 

'^'^L2(L + K.ni + K2f{L + KsfiL + K.^L + K^)' ' ^^'^^ 

Because the external momenta are "twice" linearly dependent in four dimensions, 
we can always adjust constants oq through 05 such that 

5 

aoL^ + Y,aiiL + Kif = 1, (5.4) 

i=l 

which allows us to split the hexagon integral into six pentagons. Notice that this 
depends on the dimension in which the external momenta live; dimensionally 
regulating the loop momentum does not change this. 

By using these types of integral reduction, one can arrive at expressions for 
the wanted amplitudes which contain a reduced set of integrals that are possible 
to compute, or more often, look up. 

5.1.2 Massless Integral Basis in Four Dimensions 

For four dimensional theories with massless states running in the loop, there 
is an additional identity that allows pentagons to be written in terms of lower 
point integrals [36]. This leaves a basis for all one-loop amplitudes consisting of 
scalar box integrals, triangle integrals, bubble integrals, and 10. The integrals 
are defined as (implicit ze's) 



14(^1,^2,^3) = -i/x' 



2e 



(27r)4-2. 



X 



L^L + Pi)2(L + Pi + P2)2(L + Pi + P2 + P3)2 ' 

d^-^'L 1 

(27r)4-2e L2(L + Pi)2(L + Pi + P2)2 ' 



'^^P) = -^^''j (t)ItL2(L + P)2 ' ^^•^) 

These integrals (and 1) have coefficients which are rational functions of the ex- 
ternal spinors and momenta. In addition, box and triangle integrals come out as 



^Remember that massless tadpoles vanish in dimensional regularization 
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expressions of mass dimension — 2e, divided by their respective Gram determi- 
nants which may or may not be rational. This prompts us to define the box and 
triangle integral functions F4 and where the Gram determinant is taken into 
the coefficient instead, where we will often find that it is cancelled. 

Box and triangle integrals are classified according to which of their external 
momenta are massless. For triangles there are zero-mass through three-mass (Om, 
Im, 2m, 3m) while for boxes there zero-, one-, three-, and four-mass together 
with the so-called two-mass-easy (2me) and two-mass-hard {2mh) configurations 
where the massless corners are opposite and adjacent, respectively. The integral 
functions and Gram determinants are given in appendix O 



5.1.3 Supersymmetric Decomposition 

When calculating gauge theory amplitudes, the answer depends on the number 
of massless fermions and scalars in the gauge theory because they can run in 
the loop. There are particular choices for these numbers that give dramatically 
simpler answers, namely those corresponding to supersymmetric theories. It is 
standard to use the "basis" for the particle content consisting of an = 4 
multiplet (1 gluon, 4 fermions, 3 complex scalars), an = 1 chiral multiplet (1 
fermion, 1 complex scalar), and the so-called AA = multiplet containing just a 
complex scalar. A theory with nj adjoint fermions and Ug adjoint scalars can 
thus be written as 



(AA = 4) 



n 



{M=l) + 



I - Uf + ns 



(5.6) 



If the fermions or scalars are in the fundamental representation, Uf^g should be 
substituted with nj^g/Nc where Nq is the number of colours. 

In background field gauge, it can be shown explicitly that there are diagram- 
by-diagram cancellations that remove four powers of loop momentum in the nu- 
merator for = 4 and two powers for M = 1. When those results are Passarino- 
Veltman reduced, they can only give rise to boxes in the AA = 4 case and only to 
boxes, triangles and bubbles (that is, no additional rational terms) in AA = 1. For 
J\f = there are no such cancellations, but the calculation is made simpler by the 
reduced number of Lorentz indices for a scalar running in the loop as compared 
to a gluon. 

Supersymmetric decompositions can be used in other theories as well, how- 
ever, without a formal handle on the cancellation of numerator loop momentum. 
Since the expression of the supersymmetry at tree level is the same as in gauge 
theory {e.g. SWI's, cf. section [2.3.2p and tree level results are recycled in most 
loop calculations, it is generally believed that the same regularity holds for other 
theories, in particular gravity. For arbitrary J\f it is believed that an even M 
multiplet will have Af cancellations of loop momenta in the numerator while an 
odd Af multiplet will have Af + 1 [21]. 
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5.1.4 Limits and Singularities 



To aid in both checking and construction expressions for one-loop amphtudes, 
it becomes important to look at these expressions in certain limits. The tree- 
level limiting behaviour was studied in section 12.3.31 and this is extended here to 
one-loop level. 

The multiparticle poles of Yang-Mills amplitudes can occur in the same way 
as at tree-level when a non-loop propagator goes on shell and splits the diagram 
in two, the loop being on one side or the other. We may also have a singularity 
coming from other parts of the calculation which has a more messy interpretation, 
though. This gives the multiparticle limit [30] 



E 

h=± 



+ 



p2 

_4tree(___^p/i)__4tree^_p 



p2 



(5.7) 



as — > 0, where is a function depending (in a cut-containing and non- 
factorized way) on all incoming momenta, but not on the helicities. It also de- 
pends on the particle content of the theory. 

Another type of singularity that Yang-Mills theory has, is coUinear singular- 
ities. Again, we get something similar to the tree-level case, only now we have 
to introduce 1-loop splitting amplitudes which depend on the particle content of 
the theory. We have 



h=± 



Split*^''^"(z, a^\b'"')A^-^''°^{p-^, . . .) 
+Split;^"^°°P(z, a'^", 6'''')^*'^^^(p-^ . . .) 



(5.8) 



as a ^ zP, 6 — > (1 — z)P. The one-loop splitting functions for both gluons and 
fermions are given in [33] (and to all orders in e in [32]). The pure-glue splitting 
functions are 



where 



Split i'lT 


-'°°P(a+, 


6+) = 


Spliti" 


-'°°P(a+, 


5+) = 


Split 

e 


-'°°P(a±, 




U = - 






e2 V^(l 


-z){- 



1 



1 



N N 



487r2 
crSplit*i*^'^(a+,6+) 



U- 



1 . 

1 



N N 



z(l 



L.tree i 



vr 



+ 21og2;log(l - z)- — 
b 



(5.9) 



(5.10) 
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and cr is defined in fOl) . 

The general form of infrared divergences in Yang-Mills is also known from 
general arguments. This imposes the constraint that [71,92,105] 



i=l 



+ O(e0). 



(5.11) 



This imposes some quite stringent constraints on coeffiecients of box and triangle 
functions because there needs to be complete cancellation of e.g. terms of the 
kind 



a2 



I... 



(5.12) 



where t... is any multiparticle kinematic invariant [31]. In one-loop gravity, a 
similar relation holds [81], 



_y^l-l00p 



(4vr) 



Ei<jSiilog[-s,j] 



2e 



(5.13) 



with similar consequences [48] to be explored further in chapter [71 



5.2 Unitarity Cuts 

As mentioned in the introduction, the concepts of unitarity and analyticity of the 
S-matrix came to play a large role again in the 90's. That was primarily because 
the insights of the previous section were sufficient to fill in the missing parts of the 
puzzle, the parts which (helped by experimental evidence) sent S'-matrix theory 
out in the cold. 



5.2.1 Ordinary D = 4 Unitarity Cuts 

Instead of trying to directly reconstruct the one-loop amplitude from its cut 
discontinuities — a procedure which is often ambiguous — it is possible to use all 
the insights of the previous section to come up with an ansatz for the result 
and match cuts on both sides of the equation to arrive at an answer. This was 
first used by Bern, Dixon, Dunbar and Kosower to compute the one-loop MHV 
amplitudes in = 4 Yang-Mills [33] as follows: It is known that the = 4 
amplitude only contains boxes, so we can write 

-4MHr(^">r)= E (5.14) 

aSboxes 

On both sides we can pick out the discontinuity in by replacing the loop 

propagators between ra\ and m\^ \ and between m2 — 1 and by delta functions 
in the square of the loop momenta. On the left side of (j5.14p we get an integral 
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over the product of two tree amplitudes, possibly summed over internal helicities 
in the loop, 




xA{m2 + 1, . . . , mi - 1, , {-hY^' . . . , ms, 1^ , (-/i)"^)- (5.15) 

Pictorically, this becomes 




where the dashed line denotes the cutting of the two internal lines. In Lorentz 
invariant gauges such as Feynman gauge, we would have to include a ghost loop, 
but by using the spinor-helicity formalism, which is implicitly in a light-like axial 
gauge as mentioned at the end of section 12.21 we can avoid them. 

For the two tree amplitudes to be non-zero, they must both be MHV. Then, 
if i and j are both on the left [right] side we have {h,h) = (1, —1) [(—1, 1)] while 
if they are on opposite sides h = h can take any value in the = 4 multiplet 
and we need to do an actual summation. These two distinct cases can be drawn 
as 




mi 



m2 + 1 



m2 + 1 




The summation over the = 4 multiplet, however, gives exactly the same result 
as when both i and j are on the same side, namely that the integrand can be 
written as 

l)mi){hl2) (m2(m2 + l))(^2^i) 



A 



tree.- 



((mi - l)/i)(/imi) (m2/2)(^2(m2 + 1)) 
This can again be rewritten as 



(5.16) 



^ ^tree 

2"^ 



p2 p2 

-'^rrii ,m2 mi — l,m2+l 



p2 p2 

mi — l,m2 mi ,m2+l 



(/i + (mi-l))2(/2-(m2 + l))^ 
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p2 p2 _ p2 p2 

mi,r)i2 — 1 rrii — l,m2 Tni — l,m2 — l m\,m2 

(Zi + (mi-l))2(Z2+m2)2 

p2 p2 _ p2 p2 

mi+l,m2 mi ,m2+l mi ,m2 mi +I,m2+1 



ai-mi)2(/2-(m2 + l))2 

(5.17) 



p2 p2 _ p2 p2 

mi+l,m2 — 1 mi,m2 mi,m2 — 1 mi+l,m2 



(/l - mi)2(/2 + m2)2 

plus terms that vanish when we integrate. In this way all the dependence on the 
loop momenta has been confined to two propagators in each of the four terms. 
That, however, is exactly what one would get by cutting four particular 2-mass 
easy boxes. And even better, the coefficients (apart from A^'^^^) are the respective 
Gram determinants, so this is a sum of cuts of four box functions (-F's) with the 
same coefficient. In equations. 



^mi,m2'^ ^ — ,m2 li -fmi ,m2 ) ^2 ~l~ 1) 

+F4(mi - l,Pmi,m2-l,"^2) 
+F4{mi,Pmi+l,m2,m2 + 1) 

+F4{mi,Pmi+l,rn2-l,m2)). (5.18) 



This suggests that the one-loop amplitude can be written roughly as the tree 
amplitude times the sum of all 2-mass easy box functions. The 'roughly' comes 
in through careful considerations of which choices of rrii and m2 are allowed and 
independent. The result can now be checked using other means such as collinear 
limits and it does indeed hold up. 

Staying inside Yang-Mills theory there are two extensions we can make on 
this [34] . Firstly, we can consider the = 1 or the M = theory where triangle 
and bubble integrals are introduced on the right side of (|5.14p . This complicates 
the calculation by giving integrands whose numerators depend on the cut loop 
momenta while there may be one less propagator. This then requires integral 
reduction a la Passarino-Veltman (section l5.1.ip with some simplifying conditions 
on the loop momenta. For the AA = theory there is the additional complication 
that the rational terms are not picked up by four-dimensional unitarity cuts. As 
we will see later, it can still be quite useful to know the cut-containing parts of the 
amplitude. Secondly, we can go beyond MHV amplitudes where the integrand 
does not immediately have the nice decomposition into boxes as e.g. in ()5.17p . 
Still, some AA = 4 six-point results can be obtained. 

Inspired by the recent revival of unitarity arguments, some improvements 
have been made on the technology which allows for quite efficient extraction of 
integral coefficients [58,63]. 

5.2.2 Generalized Unitarity 

The process of computing a unitarity cut where two propagators are put on-shell 
has a clear interpretation as the discontinuity across a cut in the amplitude as a 
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complex function of the kinematic variables. In the way it is applied, however, 
this interpretation of the cut is not used. Although the interpretation of the 
cutting of more than two loop propagators may not be as clear, it seems just as 
reasonable to use it to constrain an amplitude. 

The problem is that many integrals cut on more than two propagators turn 
out to be zero because there are no solutions in Minkowski space for the cut 
momenta. A triple cut is only non-zero when applied to three-mass triangles and 
three- and four-mass boxes, while a quadruple cut is only non-zero for four-mass 
boxes. Higher cuts are zero because they over-constrain the four-dimensional loop 
momentum. There are some situations where generalized unitarity does come in 
handy in Minkowski space; the NMHV M = 4, amplitude has been calculated to 
all n in this way [31,37]. 

Since the comments here apply only to amplitudes in real Minkowski space, 
they will be altered completely when we move to complex momenta below. 

5.2.3 Unitarity in D = 4 - 2e 

If we want to know a full one-loop amplitude we must do the dimensional reg- 
ularization properly. This involves taking the loop momentum to be 4 — 2e- 
dimensional rather than four-dimensional. This poses a very deep problem when 
using the spinor helicity formalism as it is tied to exactly four dimensions. One 
solution could be to use spinor helicity notation for all non-loop quantities and 
ordinary Lorentz notation for loop quantities. This comes at the expense of the 
compact form of tree amplitudes. 

When computing gluon scattering amplitudes, the easiest way to proceed is 
to use the fact that only the = component can have rational terms, and thus 
we only need the less compact amplitudes for scalar particles in the loop [26,35]. 
In this way issues of summing over polarization states are also avoided. The cut 
momentum of the scalar is written as 

L{A-2e) = L(i) + //e(_2e), (5.19) 

where L(4) is a four-dimensional momentum of mass and e(_2e) is a — 2e- 
dimensional unit vector. The integration will then split into an integral over 
a massive scalar and an integral over the mass. Using the unitarity method, this 
can be used to write the amplitude in terms of loop integrals whose rational parts 
as e — > are known. Recent versions of such methods are developed in [3,4,62]. 

5.3 MHV Rules for Loops: The BST Prescription 

When the MHV rules were conceived, the perspective for their use at loop level 
was rather dim. It was argued that the loops would necessarily contain states of 
conformal supergravity, a theory which is non-unitary. Not everyone was deterred 
by this argument, at least not Brandhubcr, Spencc and Travaglini (BST) who 
decided to recalculate the one-loop A'" = 4 MHV amplitude using MHV rules [57] . 
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Their prescription consists of drawing all diagrams with two MHV vertices 
and two internal lines to form a loop, 




The momenta Li,L2 of these internal lines are off-shell, so on-shell continuations 
of them are required for use in the MHV vertices. This is performed by taking 
the integration measure 



and changing variables to li^2 and 2:1^2 such that li^2 are lightlike and 

Li,2 = h,2 + zi,2V (5.21) 

where rj is an arbitrary lightlike vector. As shown in [55], the combination zi +Z2 
can be integrated out by standard (albeit involved) methods of complex analysis, 
while the combination zi — Z2 can be replaced by = (P + z\ri — 22^)^- This 
yields the measure 

2md%d%5M (;2) ^f^^^p^ ^(4) (/, _ + P,). (5.22) 



where the on the delta functions impose that li and I2 must have positive 
energy. 

For = 4 matter running in the loop, BST found that the product of the 
two MHV vertices (summed over the multiplet) became — in a way similar to the 
unitarity example above — a sum over terms of the kind 

t,ee jPz + mifjP, - m2f - PKPz + mi - m2? 

2{h + rmY{l2+m2Y • ^ ' 

At this point the loop momenta are described as 4-momenta rater than spinors, 
and that permits a dimensional regularization by taking the number of dimensions 
to be 4 — 2e as usual. 

With sufficient hard work, the integration of (j5.23p with the measure (j5.22p 
can be performed and terms (four of each) which have the same numerator in 
(I5.23P combine to form two-mass easy box functions. The proof relies in part on 
an identity between nine dilog functions and in this way gives a different form of 
the two-mass easy function derived earlier by Duplancic and Nizic [82] but not 
widely known. 
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There are generalizations to the cases with J\f = 1 [13, 124] and scalar matter 
[11] in the loop which require triangle and bubble functions. In the scalar case 
there are additional rational pieces which are not predicted by this version of 
the loop MHV rules because they contain cuts almost by construction. These 
rational terms, however, appear to be generated in the light-cone reformulation 
of the MHV rules with a particular non-dimensional regularization [56] . 

In spite of the successful application to several one-loop calculations of MHV 
(and similar [6,10]) amplitudes, the EST prescription appears to be unpractical 
for higher numbers of negative helicity gluons than two. Although the formal 
correctness of the prescription has been proven using the Feynman Tree Theo- 
rem [55] for any number of negative helicity gluons, there have been no concrete 
calculations, mainly because the method has been superceded by others men- 
tioned in this chapter. 

The prescription also inherits a formal problem from the tree-level MHV rules, 
namely the one mentioned in section 13.2.21 that the vertices do not obey mo- 
mentum conservation. From that perspective, it may be seen as odd that the 
right results appear, as any (momentum conserving) rewriting of the vertices will 
change the results to something wrong. Resolution of this problem comes either 
from formulating the loop-level MHV rules in a light-cone form or by noting that 
any cutting of internal lines (as in the Feynman Tree Theorem mentioned in the 
introduction) reduces the amplitude to integrals over on-shell tree amplitudes 
where the conservation-of-momentum problem was resolved by viewing the rules 
as coming from recursion relations in complexified Minkowski space (chapter H]) . 



5.4 Quadruple Cuts 

The unitarity method of section [5^2] employed cuts on two internal states, partly 
because the cut conditions are guaranteed to have solutions for which the internal 
momenta are real in Minkowski space. Clearly, there is nothing wrong with 
applying more cuts in principle, but the cut conditions often turn out not to allow 
any real solutions, thereby setting the cut to zero and rendering its information 
useless. 

While working with twistor methods, first described in 2-1-2 dimensions where 
the spinors are real and independent, Britto, Cachazo and Feng realized [59] that 
working in that signature (or, more generally, in complexified Minkowski space) 
would allow simultaneous solutions to the constraints of four cuts. Since the 
loop momentum is four dimensional this would fix it completely, reducing the 
problem of finding the coefficient of a box function to a purely algebraic one. 
This is particularly useful for calculating = 4 amplitudes which only contain 
box integrals. 
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5.4.1 Basic idea in jV = 4 



Consider first a one-loop amplitude and choose four internal lines to take on-shell. 
There is, of course, a one-to-one correspondence between the boxes and the sets 
of four internal lines, so there is no double counting issue. The chosen box with 
coefficient c is 



^^L 1 

c 



4 [L2 + ie][{L + Kif + ie][{L + Ki + + ie\[{L - K^f + ie] 

(5.24) 



which is cut to 



j d''L5^+\L^]8^+\{L + K^f]5^+\{L + K2f]8^+\{L-K^f]. (5.25) 

soln. 

The constraints have two solutions which give rise to the same Jacobian for the 
integral. The Jacobian is inversely proportional to the Gram determinant but 
will otherwise not concern us. The result of the cutting is now 

2c:i&c{Ki,K2,Ki). (5.26) 

We can do the same considerations for each Feynman diagram in the calculation 
of the one-loop amplitude. Adding up contributions from all diagrams that cut 
our four chosen internal lines will give 



/ 



d^L J2 A{...,-L,L + Ki)A{...,-L-Ki,L + Ki + K2) 

soln., states 

x^(. ..,-L-Ki-K2,L- K4)A{. ..,-L + K4, L) 

x,5W [L^] <5W [(L + Krf] [{L + i^s)'] [(L - K^f], (5.27) 



where the summation over states refers to the possibility that different particles 
may have those momenta, and that we must sum over all possibilities. In the 
same way as above, this gives 



f AAABAcADyac{Ki,K2,K4), 

soln., states 



(5.28) 



or, 



c=\ ^aAbAcAd- (5.29) 



2 

soln. ,states 



In other words, the coefficient is given by products of four on-shell tree amplitudes 
in four dimensions. The quadruple cut can be depicted as a box with four corner 
tree amplitudes connected by the (implicitly cut) loop propagators. 



61 




In concrete examples the internal momenta have helicity labels if there is only 
one helicity configuration which gives non-vanishing contributions. 

5.4.2 A Simple Example 

As a first example, we may try to calculate the one-loop ( amplitude. 

This has only one quadruple cut, namely 



b- 





^3 






^ + 




± 
















± 










h 





or d+ 



There are two solutions for the internal momenta l\ through ^4. The first is 



,3 = |*)||^[c|, '4 = M)||[c|. (5.30) 

with the consequence that MHV three-point amplitudes at corners B and D are 
zero while googly-MHV three-point amplitudes at corners A and C are zero. 
Thus, the only non- vanishing assignment of internal states is 
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giving the contributioijl 



{ah 



d+ 



M-h)? {ki-h)f [i-h 



1 

2jl2{-h)){{-h)a) [i-hMbh] {{-h)c){ch) [<i/i][/i(-/4)] 
1 {ahf hk? ihkf [hdf 



2{ah){hl2) mibh] {hc){ch) [hhWhd] 
l{daf{dcf{bcf {abf [acf {bdf [cdf 
2 {ah){lib)[ab][bc]{bc){bd)[cli][h 
l{daf{dc)^{bcf {ab)'^ [ac]3 (bd)'^ [cd\^ 
2 



-\stA 



{bd)^ab)^ 

tree 



{ad){db) [ab][bc] {bc){bd) [ca][ad] 



(5.31) 



where s and t are the normal Mandelstam variables. The other solution for the 
internal momenta comes about by 'flipping' |-) <-> |-] in (|5.3U|) and appropriately 
reversing the conclusions about the helicity assignments to get 




Notice how the symbols © and © are used for three-point MHV and googly-MHV 
corner amplitudes, respectively. Since A^^^^ is invariant under flipping we get the 
same result as before. The one-loop amplitude then becomes 

^l-loop ^ _st^treej^(^^ ^) ^ 2A''^''Fi{a, b, c). (5.32) 



Remember that |( — fc)) = and |(— fc)] 



|fe]. 
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In practice, a calculation will be a good deal more involved as it may include more 
complicated corner amplitudes, especially non-MHV amplitudes, and summing 
over the full N = A multiplet possibly running in the loop. Such an example 
(pieces of which were calculated in the original article on quadruple cuts [59] ) is 
presented in the next chapter. 

To aid our calculations in the remainder of this thesis, the solutions of the 
internal constraints have been assembled in Appendix [Bj 



5.4.3 Box CoefRcients in Other Theories 

The quadruple cut method extends smoothly to other theories, but without be- 
ing "complete" in the sense that the entire amplitude can be found from that 
method. What it does calculate is the contribution to the amplitude from box 
functions, since the quadruple cut singles out the boxes. When calculating a gen- 
eral Yang-Mills one-loop amplitude, the contributions to boxes from the M = 1 
chiral multiplet running in the loop can be calculated just by adjusting the in- 
ternal states. As is done in chapter [HI the method may also be used where some 
(massive) external states do not run in the loop. Another extension is to gravity 
theories such as = 8 supergravity, where it is possible in general to use the 
KLT relations between tree amplitudes to construct KLT-like relations between 
the box coefficients. The use of quadruple cuts in gravity will be considered in 
more detail in chapter [71 



5.4.4 GeneraHzed Unitarity in Complexified Minkowski Space 

As briefly mentioned in section [5.2.21 taking quadruple cuts was not a genuine new 
idea in itself, rather it was the use of internal momenta in complexified Minkowski 
space. Thus, the method also allows for wider use of the triple cut. This leads 
to different methods for obtaining the cut-containing parts of amplitudes [51,85, 
88,118]. 

Also, both triple and quadruple cuts can be done with massive particles or, 
equivalently (c/. section [5.2.3p with massless particles exactly in 4 — 2e dimensions 
to yield also the rational parts of amplitudes [54, 114] 

The idea of generalized unitarity can be extended to higher loop orders, where 
cutting a large portion of the propagators can reveal information as to the struc- 
ture of both the momentum-dependent numerator and any other coefficients, and 
can likewise provide useful checks of higher loop unitarity cut calculation where a 
more conventional generalized unitarity approach is taken. An amusing case [64] , 
which unfortunately does not seem to generalize, is = 4 two- loop integral co- 
efficients up to and including six external gluons where completely cut integrals 
such as 
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have an additional propagator-like singularity coming from the Jacobian of the 
right box. This singularity can again be cut providing a total of eight constraints 
which fix the integration completely and gives the coefficient as a product of tree 
amplitudes. This, however, requires that there is a box which is two-mass hard 
or one-mass; otherwise the Jacobian is not propagator-like. 

Recently, Cachazo and Skinner has shed some additional light on the use of 
cutting procedures at higher loops [66] . 



5.5 Recursion at Loop Level 

Up till now, this chapter has been concerned with deriving one-loop amplitudes 
from tree amplitudes. In that perspective it should come as no surprise that 
recursion can also play a role here. Indeed it does play a role in determining 
coefficients of the integral functions [29,51], although it does not seem to be as 
strong a tool as at tree level. 

The problem of determining integral coefficients is, however, not our main 
objective in this section. Rather, it is the calculation of the rational parts of 
non-supersymmetric one-loop amplitudes directly in four dimensions without ap- 
pealing to D-dimensional (generalized) unitarity. The methods described here 
were developed by Bern, Dixon and Kosower, first for amplitudes with no cut- 
containing parts [38,39] and later for amplitudes with cut-containing parts [40]. 
Since this section is more intended as background material for chapter [H] than a 
review of advanced topics, the treatment will be adapted to that case and other 
complications will only be mentioned briefly. 

As in the tree case, on-shell recursion starts out by deforming two external 
momenta with a complex parameter z and writing the 2 = result as an integral 
along a closed contour enclosing the origin, 

^i-'°°P = ^ / — ii-i°°P(z). (5.33) 

^'^'^ Jaround ^ 

If ^^~^°°P(z) — > as z ^ oo, we can write this formally as a sum over integrations 
around all singularities of „4,^~^°°p(2:), which may be both poles and branch cuts, 



2m 

singularities i 



E 



around i 

ResiA^^-i°°P(z) f dz 



1 • 

poles I 



[ — DiscB^^-i°°P(z) (5.34) 
Jb z 
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where i? is a branch cut and Disc^ is the discontinuity across it. In the first term, 
the residues are given by the appropriate multiparticle or colHnear factorization 
and depend on lower-point one-loop amplitudes or splitting functions. Since we 
assume that the cut-containing parts of the amplitude can be calculated by other 
methods described in this chapter, it makes sense to split all one-loop amplitudes 
and one-loop splitting functions into their cut-containing and rational parts. 
Then it can be argued that cut-containing parts can be deduced from the cut- 
containing parts of lower-point one-loop amplitudes and splitting functions, and 
that rational parts can be deduced from the rational parts of them. 

This removes the complication of cuts in the amplitude, but introduces a new 
one. When we make the split of an amplitude into cut-containing and rational 
parts 

^i-ioop = c + R, (5.35) 

we cannot trust that C and R individually only have physical {i.e. multiparticle 
and collinear) singularities. In fact, the cut containing parts will often contain 
terms like 

(5.36) 

which do not correspond to physical singularities of ^^~'°°p. To handle this, we 
add and subtract a rational cut- completion term CR, 

^i-ioop = ^c + CR) + {R- CR) (5.37) 

which has (minus) the same unphysical singularities as C. The singularities of 
R — CR are poles (single or double) in kinematic invariants, so we can do the 
usual on-shell recursion trick (again, assuming R{z) — CR{z) — > as z ^ oo), 

1 ■ 

poles I 

where the poles are those that arise from known sing ularities of ^i-'°°p. The 
terms coming from R and CR are known as the direct recursive DR and overlap 
O terms, respectively. The latter can be computed directly from CR{z) while 
the former is more intricate. 

If the one- loop splitting functions had only single poles, we could interpret 
them — in the same way as tree splitting functions — as recursion diagrams of the 
type 



Rational terms proportional to tt are included with the cut-containing part. 
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Drawing it like that, we can now use the rational parts of the splitting functions 
(|5.9p for a gluon it the loop, 

Spliti-'°°P'^(x,a+,6-) =0, 
SplitV''°°P'^(x,a+,6+) = -^7^(r3^i«^ 



3 ' '{ab)^' 

Spliti-'°-'^(x,a+,6+) = ""^ (5.39) 

3 {ab} 

to say something about the contribution of such a diagram. Firstly, if a and b 
have opposite helicities, there is no contribution. This is also the case if all three 
helicities sitting on the three-point "vertex" and a had its antiholomorphic spinor 
shifted, since the requirement of d and b being collinear is [db] = 0; conversely, if 
a had its holomorphic spinor shifted the all-minus splitting function would give 
zero. 

If the holomorphic spinor of a is shifted and we have to deal with an all-plus 
splitting amplitude, we immediately see that it has a double pole in z. This poses 
a serious problem because finding the residue requires us to seperate the pure 
double pole from any single pole which may be hiding underneath it. However, 
Bern, Dixon and Kosower found a way to circumvent this by introducing a slightly 
different three-point vertex which seems to do the job in certain cases. Lastly, 
the complex factorization of the last splitting amplitude above is quite murky, 
and it should generally be avoided [39]. 

Keeping these matters in mind, we can obtain the direct recursive terms as 



DR= [ ^2 + ^2 



i,h=± 

\tree( T3h\ Atrccfry—h 



where i runs over all channels affected by the shift, R are the rational parts of 
lower-point amplitudes or properly defined three-point loop vertices and J-ji is 
the rational part of the factorization function. Summing up, the rational part of 
the amplitude we are after is 

R = DR + CR-0. (5.41) 
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A last complication is that of the z — > oo behaviour of the amplitudes. It is not 
always possible to find a shift which is both zero as z — > cxd and avoids problematic 
splitting functions. However, it has been shown [19] that all problems can be 
circumvented with imaginative combinations of shifts. This means that any one- 
loop amplitude in massless QCD is in principle obtainable by the combination of 
unitarity methods for the cut-constructible parts and recursion for the rational 
part. A computer program under the name of BlackHat which performs such 
calculations is under development [18]. In chapter [8] we will turn to a complete 
calculation of a particular one-loop amplitude in QCD with an effective coupling 
to the Higgs particle. 
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Chapter 6 

One-Loop Amplitudes in AT = 4 
Super- Yang— Mills 

In this chapter we take a closer look at the calculation of one-loop amplitudes in 
A/" = 4 Yang-Mills. As described in section EH these amplitudes can be written 
as linear combinations of scalar Feynman integrals, where the coefficients are 
essentially products of four tree amplitudes. Here we will calculate some concrete 
one-loop amplitudes as a demonstration of the technique. The first calculation, of 
the one-loop MHV amplitude, is primarily included for demonstrational purposes. 
The part of the calculation of the one-loop NMHV amplitude where scalars and 
fermions in external states are taken into account is based on [127]. 

6.1 MHV Constructibility 

Although the quadruple cut method allows for the efficient numerical evaluation 
of one- loop amplitudes from tree amplitudes, there are many analytic simplifica- 
tions to exploit. The most important of these is the combination with the simple 
all-n expression for MHV tree amplitudes. For certain configurations of external 
helicities, and for some box coefficients, the required tree amplitudes are all ei- 
ther MHV or googly-MHV, leading to easily obtainable all-n expressions. A box 
coefficient with only MHV and googly-MHV corners is called MHV constructible. 

The simplest class of MHV constructible box coefficients are those of the MHV 
one-loop amplitude. They were first computed by the unitarity method of section 
I5.2l by Bern, Dixon and Kosower [33], but we will review the calculation from the 
quadruple cut perspective here for the sake of familiarity with the method. 

As in any quadruple cut calculation, we must start by determining which box 
functions contribute. This can often be reduced to a matter of counting pluses 
and minuses: The MHV amplitude has two external minuses. Inside the cut 
box there are four internal lines with one minus each. Thus, the four tree-level 
corner amplitudes must have a total of six minuses. The only way these can be 
distributed without giving zero amplitudes is to have two MHV corners and two 
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googly-MHV three-point comers, and the latter must be at opposite corners. In 
other words, we can only have 2-mass easy boxes where the massless corners are 
googly. Of course, this argument only holds for gluons in the internal and external 
states, but the conclusions hold in general because all non-zero amplitudes are 
related to the gluons-only amplitude with the same sum of helicities. 




Rather than now classifying the positions of the two negative helicity external 
legs, we can start out by using the simple form of the MHV amplitude to write 
down part of the coefficients, namely the part that depends on the denominator 

of dzi]), 

1 1 1 

2 {{-h)d +i){d+i...b-i){b- ih){i2{-h)) [{-hmhhWH-h)] 

1 1 

{{-h)h + l){b+l...d-l){d- l/4)(/4(-/3)) [{-l^)d][dh][h{-U)\ 

_ 1 {dd + 1) (fe - lb) {hh + 1) {d- Id) 
~ 2(- • •) [dh] {hd +l){b- 1/2) m m {kb +l){d- Ik) [kd\ 



{hh)[l2h]{hk)[kh]' ' ■ ' 

Using the solutions of the internal constraints ()B.3p , we can write this as 

1 {db) {bd) (bd) (db) {bd)^ 
2{- ■ ■) [dAb) {dAb] [bCd) {bCd] {db){dACd){bd){bCAb) 

_ {bcd\{dAb] f {bd) y 

2(---) \{bCd]{dAb]) ■ ^ ' 

Notice that 

{bCd]{dAb] 
= -{bC{A + b + C)Ab] 
= -A^{bCb] - C'^{bAb\ - {bAb]{bCb] 
= A\C^ - (6 + Cf) + C\A^ -{A + 6)2) 

-{A^ -{A + bf){C^ -{C + bf) 
= -{A + bf{C + bf + A^C^ (6.3) 
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is twice the Gram determinant. Thus, the denominator consists of the Gram 
determinant, the conventional (• • •) spinor product, and a remaining term to the 
fourth power. 

The term coming from the numerator depends on the position of the two 
external negative helicity legs. These can be placed in four distinct ways, but we 
will content ourselves with the calculation of two of them. The first case is when 
the two negative helicity legs are on the same MHV corner, say corner A, 




If we call them i and j, respectively, the numerator becomes 



,...4f[bCd){bd){bCd]Y 

V'J) 



(6.4) 



We see that the paranthesis cancels with the one in (16. 2p and leaves only the 
Gram determinant times the tree amplitude. 

We now turn to the case where i and j are on corners A and C, respectively. 
In this case, the whole N = A multiplet can run in the loop which we draw as an 
extra internal line. 
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The numerator contribution from the gluons can easily be seen to be 

{il2)'[{-l2W{jh)'[{-k)d]' + {{-h)i)^[bk]'{{-h)j)'[dh]\ (6.5) 

The contribution of the fermions can be computed by noting that the amphtude 
A^'^^{g-,f-,f+, . . .) is simply related to A^'^^ {g- , , . . .), namely by the ex- 
change of one of the numerator factors with another. This yields 



-4[{ii2y{i{-h)) [{-hn'm {juyiji-h)) a-rn'm 

+ (/2.)((-/l)^)3[6(-/2)][W3]='(/4j)((-/3)jfK-/4)][<i/l]3), (6.6) 

where the prefactor takes into account that there are four species and that they 
must have fermionic sign. The scalar contribution becomes 

6(i/2)'(i(-/i))'[(-Z2)6]2[/36]'0-^4)'0-(-^3))'[(-/4)^^]'[/i^^]'. (6.7) 

When we add these three expressions together, we see a common situation in 
= 4 SYM: the particle multiplicities make sure that the result is simplified, in 
this case by making it writable as a fourth power of some term: 

\ii2)[{-i2mjk)[{-k)d\ - {{-i,)i)[bk]{{-k)j)[dh]"' 

m[hmi)[hd]-{Ui)[hd]{jk)[hbY^ 

{ij)^{{hh)m[hd]y 

^ _ / \ni ;n ) ( nn w n ai /7 1 \^ 



f [bCd){bd){bAd] \ 
\ ) 



In words, the particle states of = 4 running in the loop conspire to make the 
numerator factor the same as before. It is not difficult to imagine (and true) that 
all four positionings of the negative helicity legs i and j give the same numerator 
factor. 

We can now write down the complete one-loop MHV amplitude in = 4 
SYM. It is 

n min(n,i— 2+n) 

<"Hr = cr<lvE E F^{hn+i,-i.j)- (6.9) 

i=l j=i+2 

This was the result for gluons in the external states, but in fact, it holds 
for any MHV choice of external particles. Remember from section 12.3.21 that 
SWI's fix all tree-level MHV amplitudes up to a constant; since the one-loop 
result should also respect the full = 4 supersymmetry, the same is the case 
here. Only the common factor is corrected from its tree-level value to its one-loop 
value. 
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6.2 NMHV Amplitudes 



The next natural step would be to compute all NMHV amplitudes. The purely 
gluonic case was worked out by the unitarity methods of section 15.21 by Bern, 
Dixon and Kosower [37] prior to the publication of the quadruple cut method. 
With this method, its inventors computed the MHV constructible parts of the 
NMHV amplitude [59] . The case with other external particles was worked out by 
Bidder, Perkins and the author [127], but since the latter calculation naturally 
incorporates the pure gluon case, both cases will be reviewed here. 

For most (but not all) of the box coefficients of the NMHV amplitude, we 
can again exploit MHV constructibility. Doing the same analysis of the "number 
of minuses" as above, we see that there are a total of seven minuses, which can 
be distributed among the four corner amplitudes in two ways: Either there is 
one googly-MHV three point corner and three MHV corners, or there are two 
googly-MHV three point corners, one MHV corner, and one NMHV corner. The 
latter solution is not MHV constructible and will be postponed until we have 
considered the former which is. The diagram for the calculation and the solution 
for the internal spinors are 




\d){dCBA\ 

{dACd) 
\BCd){dA\ 

{dACd) 
\BAd){dC\ 

{dACd) 
\d){dABC\ 

{dACd) 



(6.10) 



6.2.1 3 Mass Boxes, Gluons Only 

As in the MHV case, we start by getting the denominator factors out of the way 
since they are independent of the external helicities. To fix the ambiguity in scale 
of the spinors, we assume the convention 

|/i> = |'i), \h) = \BCd), \h) = \BAd), \k) = \d). (6.11) 

The numerators now become 

1 1 

2 {{-l^)Ai){Ai ■ ■ ■ A^i){A^ih){l2{-h)) 

1 

{{-h)Bi){Bi ■ ■ ■ B_,){B_ik){h{-l2)) 
1 

{{-h)Ci){Ci ■ ■ ■ C.i){C-ik){h{-h)) 
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X 



(6.12) 



[{-l4)d\[dh]M-h)] 
1 (d^i)(A_ii?i)(i3_iCi)(C-id) 
2(- • •) {Aih)[hU]{kC.i)[dU]{hk){l2h)[hd\ 
1 

1 -{dACd)^ 
2(- • •) {dCBACBAd) 

[dCh){hAd\{A.ih){hBi){hh){B^ih){hCi) 

2(---) (tiCS^d] [dCBAd){dCBAd]{dCBBAd) 
1 

^ {A_iBCd){dCBBi){B^iBAd){dABCi) 
jdCBAd] 1 
2 (dC5^d]4 

^ -(^-ii3i)(i?_iCi) 

(• • ■)P^{dCBA^i){dCBBi){dABB^i){dABCi) ' 



By manipulations of the same kind as (16. 3p it can be shown that {dCBAd]/2 is 
the Gram determinant. We see the same kind of structure as for MHV: the Gram 
determinant, something to the power 4 (which we hope disappears) and a factor 
involving (• • •) which looks natural as a box coefficient. 

We can now turn to the numerators. The easiest way to proceed is to consider 
the gluon case first and then use NMHV SWT's and the known facts about the 
MHV amplitude for other-than-gluons to relate them to the gluon amplitudes. 
There are seven distinct non-zero configurations of the three negative helicity 
gluons and we will show the calculations for three of them. 

The first case is shown in the box diagram below 




The numerator is given by 
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= - {mim2) {hU)[hli]{hl2?j 
= {{mim2){hd)[dAh)f 

= i{mim2){dABd))'^{dCBAdf. (6.17) 

We should be happy to see the factor {dCBAd]^ appearing as it cancels the 
corresponding (unphysical) factor in the denominator. The second placement of 
the negative helicity gluons we will explore is 




A d+ 



which has the numerator 

[{mim2){{-l2)m^){{-U)h)[{-k)d])^ 

= (^{mim2) ihms) ikhd]^ 

= i-{mim2){l2m3){l3Cd])^ 
= {{mim2){dCBd)f{dABCdf 

= i{mim2){dCBd))^{dCBAd]\ (6.18) 

As before, we see the wanted factor to the power four. The last case is the one 
where the three negative helicity gluons are situated on the three massive corners, 




This configuration is rather special, as the whole N" = 4 multiplet can run in the 
loop and we have to sum over all the particle species. The two gluon contributions 
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can be constructed as for the cases above, 

4 



((mi/2) {m2h) {m3k) [di-h)]) + {{mii-h)) K (-/s)) (mg (-^3)) [dh]) ' 



+ 



(^{miBCd){m2BAd){m3d)^ + - {mid){m2BCd){m2,BAd)^ 

jdCBAd]^ 
^ {dACd)* ' 



(6.19) 



As in the one-loop MHV calculation above, the fermion and scalar contributions 
follow from the gluon ones by taking products of the two parantheses above with 
powers summing to four and by adjusting the prefactors. This exactly turns the 
final result into a complete power of four: 

'{mil2){m2h){m3U)[d{-U)] + (mi(-/i))(m2(-/2))(m3(-;3))[c^^i])^ 

(^{miBCd){m2BAd){m3d) - {mid){m2BCd){m3BAd)y^^^^ 

(mami) {dBCd) {m2BAd) 

4 {dCBAdf 



+ 



{dACd) 



{mid){{m3BCd){m2BAd) - {m2BCd){m3BAd)^ 

= - {m3mi){m2BAd) + Kl{mid){m2m3)y {dCBAd]"^ 
= (^{mim2){m3BAd) + {m3m2){miBCd)y {dCBAdf. (6.20) 

Again, we see the same picture emerging as before. 

When doing all negative helicity gluon configurations, we finally obtain the 
result 



c 



3m^^- ) = -{A_iBi){B^iCi)nt 



{■ ■ ■)P^{dCBA_i){dCBBi){dABB_i){dABCi) 
where k denotes the configuration and 

TiBBd = {mim2){dACd), 
Ti-ABd = {inid){m2BCd), 
Ti-BCd = {m2d){miBAd), 
Ti-ACd = {mid){m2d)K%, 
T-Laab = {mim2){m3BCd), 
Ubcc = {m2m3){miBAd), 
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'Habb = {m2ms){miBCd), 
'Hbbc = {mim2){m^BAd), 
Ti-AAC = {mim2){mzd)K%, 
Ti-ACC = {m2m3){mid)Kl, 

T-Cabc = {mim2){m-iBAd) + {m3m2){miBCd). (6.22) 



6.2.2 NMHV SWIs 



Before proceeding to the equivalent calculations with other particle content than 
gluons, we should consider some of the constraints on those amplitudes coming 
from supersymmetry. In section 12.3.21 we looked at the effective supersymme- 
try of tree amplitudes, supersymmetry algebras and the supersymmetric Ward 
identities that follow. Since we are now working with a genuinely = 4 su- 
persymmetric theory we can use these considerations for loop amplitudes also. 
Even better, since the box functions are all independent, = 4 SWIs for the full 
amplitude imply that they must also be satisfied coefficient by coefficient. Here 
we will derive SWIs for amplitudes A where all states are implicit. 

To derive NMHV SWIs relating pure glue amplitudes to amplitudes with a 
pair of fermions, we start from A{gi ,g~ ,gj , f^^) and apply Qail, G) to find 

= (l<z)^(/-i,5-,5-,/i) + {iq)A{g^J-^,gJjXk) 

+ {jq)A9l,9^,fajJa,k) - {kq)A{g^,gr,g-,gl). (6.23) 

This relation holds irrespective of the number of supersymmetries since we are 
only using one. In a notation closer in spirit to what we will use here, the SWI is 

Q={lq)A(l-^'\r,r,k+^l^) + {iq)A{l',r^/^r,k+^'^) 

+ {jq)A{l~X.r"\k+^'^)-{kq)A{l-X.rX)- (6.24) 

We can extend this to scalars in external states by starting from the amplitude 
A{9i ,gi Jh^pSab,k) and applying again Qa{q,6)- Using (|2.49p we get 

= {lq)Aif^^^,g:rJ^ .,Sab,k) + {iq)A{g^,f~-, f^j,Sab,k) 

1 9i 1 ^ab.j 5 ^ab,k 

)-{kq)A{g^,g;,f^^^,f^^,), (6.25) 

which relates amplitudes with both gluons, fermions and scalars. For this SWI we 
used that an amplitude involving A{. . . , Sab, fa) vanishes, something which is not 
immediately obvious. This follows roughly by noting that s^fe ~ Qafil^ ~ Qbfa 
which can be interpreted in the way that Sab contains a f^, and this sets the 
amplitude to zero by fermion conservation of helicity. These arguments can be 
made more exact by considering a superfield formulation of the theory such as 
Nair's [115] which plays a major role in some twistor studies. 
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To obtain SWIs involving amplitudes with no negative helicity gluons we can 
start from A{g^ , f^^, f^j, Sab,k, fc,i) and apply Qa{q, 9) to get 

= /c7j> /foj' ^oh,k, f^i) + {iq)A{g^ ,Sac,i, f^j, Sab,k, f^i) 

+{jQ)A9i,f;,, sab,j,sab,k, ft,i) - {kq)Aig^,f-,, fb^J^k^ fti)- (6-26) 

It is possible to deduce quite a few of these NMHV SWI's and we shall not 
attempt to derive all of them here. The general procedure for generating them is 
this: Start with an amplitude which has three and a half units of negative helicity 
(N^/^MHV) and where the SUSY indices b, c and d are each maximally used once 
on /~ and Sab states, and where the indices a, 6, c and d are used maximally once 
on /+ and Sab states (and a must be used once). Then, assuming the four SUSY 
indices to be different, act with Qa{Q,S) to get the SWI. Resulting amplitudes 
involving a twice in and Sab states are zero for the reason discussed above. 

As should be clear from the discussion of SWTs, they are not concerned with 
colour ordering, or numbering in general. 

6.2.3 3 Mass Boxes with Two Fermions 

We now want to calculate the three mass box function coefficients when we have 
two negative helicity gluons and a pair of fermions. We first review some general 
arguments and then give the full results without detailed calculation. 

The easiest way of classifying the configurations is to view it as a process of 
moving a half unit of negative helicity from one of the m^'s to a leg q situated on 
one of the four corners. Take as an example the configuration q, mi G A, 771-2 G B 
and 7713 = d and move a half unit of negative helicity from mi on yl to g which is 
also on A. Since this only involves corner A the ratio of the coefficients must be 
the ratio of the A corner amplitudes. 




c3-({mr},{m2},{},<i-) 

{hq) 
{hmi) 
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{mid) 

= (6.27) 

nABd{mi,'m2) 

If the —1/2 helicity came from d we could track the effect on each corner and 
derive the relative factor in a similar way. However, we could also use the SWI 

{qv)c'^{{m^},{m,},{},d-) 

= (mi7?)c3'"({gl/2^ ^-l/2|^ 

+ (m2?7>c3-( V/2, ^-}, {m-'/'}, {}, d~) 

+ {dr^)c'"'i{q'/\m^},{m^},{},d-'/^), (6.28) 
and set rj = m2 to obtain 

c3"^({gV^m^},{m,},{},d-V^) 
c3-({m^},{m2},{},d-) 
_ {qm2) _ (mim2) 'HABdiq,m2) 
{dm2) {dm2) 'HABd{mi,m2) 
_ (mig) 
(mid) 

^ 'HAAB{mi,q,m2) ^g) 
'HABd{mi,m2) 

The case chosen is, of course, deceptively easy since c'^'^{{q^^'^ ,m^}, ^^'^}, {},d~ 
as can be seen by trying (and failing) to assign internal helicities in the quadru- 
ple cut. Taking the more involved road, it can also be found by using the SWI 
above with ij = mi. 

In most cases there are no large advantages of using the one method (tracking 
the half unit of helicity through the vertices) over the other (SWI) but in some 
cases the former should definitely be avoided. Take as an example the gluon 
configuration mi £ A, m2,m3 £ B and move —1/2 helicity from 7713 to g G C. 
When analysing the possible internal helicity assignments we find that there are 
in fact two different ways it can be done. 




A m^ 
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Rather than calculating the result from this analysis, we can use that (including 
a fermionic minus) 



B -1/2 




A 




A m 



„3m / 



c3'"({mj"}, {m^, m^}, {}, d+) 
U-h)) [dh] {{-k)q) 
ihmi) [di-h)] {i-h)k) 
ihAdlikq) 



{l2mi)[dCh) 
{dABq) 
~ ^ {dCBmi) 
_ 'HBBcim2,m3,q) 

Ti-ABBimi, 1712, 1113) 

and use this in the corresponding SWI with rj = 777-2, 



(6.30) 



„3m 



{{m^ },{m2 ,7773 },{g^/^},d+) 
c3™({7775"}, {777^, 777;^}, {}, d+) 



{qm2 



+ 



(77717772) H B BC {1^2 , , q) 



(77737722) (77737712) TiABBimi, 1112,1713) 
{qm2){miBCd) + {mim2){qBAd) 
(17131172) {mi BCd) 
Ti- ABC (mi, 1172, q) 



T-iABBimi, 1712, 1713) 



(6.31) 



We see the continuation of the picture that moving —1/2 helicity from one corner 
to another produces a ratio whose fourth power is the ratio for moving —2 helicity 
between the same legs. This holds not only when the external helicities constrain 
us to a single internal state in the loop, but holds when more internal states open 
up. The latter is a consequence of the M = 4 supersymmetry. 

This becomes a helpful realization when we want to move helicity around in 
the ABC gluon configuration. No matter how we do it, we will need to sum over 
at least eight states of the multiplet. This can of course be done (and was indeed 
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done in [127]) but it is simpler to note e.g. that 



{qv)'HABc{mi,m2,m3) = {miri)TiABc{q,m2,m3) 

- {m^rtjUAAC (mi ,q,mz) 

-{m3r])HAAB{mi,q,m2), (6.32) 

and interpret this to be 'proportional to' the corresponding SWI 

{qv)c^"'{{m:[}, {m^}, {m^}, d+) 
= (ml^?)c3™({m^^/^ gV2}^ {^-}^ {^-}^ ^+) 

+ (m2r/)c3™({mr, {m^'^^}, {m^}, d+) 

+{m3v)c'"'{{m^,q'/'}, {m^}, {m.'^'}, d+), (6.33) 



in the sense that 



c 



c3"*({mi },{m2 },{m^ },d+) Habc (mi, 1712,1713)' 



(6.34) 



etcetera. Using the methods described in this section, one can derive all box 
coefficients involving two gluinos. They are reported in the tables below. The top 
line of each table states the gluon configuration we start from and the associated 
H factor. Below, the first column states where the —1/2 unit is moved to (the 
position of q) and the second column states which of the negative helicity gluons 
it came from. The third column then states which factor the TL in the top line is 
to be replaced by to obtain the correct box coefficient. In the first three cases we 
have used identities such as —Ti.ABd{^i-i''^2) = Ti-AAB{q-,'mi,m2) to abbreviate 
results. 



AAB 


'HAAB{mi,m2,m3) 




AAC 


Ti- AAC {rni, 1712, 1713) 


A 


nil 


'HAAB{q,m2,m3) 




A 


nil 


T~iAAc{q, m2,ni3) 


Id 




'HAAB{mi,q, 1713) 




Id 


ni2 


'HAAc{mi,q, ma) 




1713 









ni3 





B 


nil 


'HABBim2,q,m3) 




B 


nil 


-'HABc{q,ni2,ni3) 




ni2 


'HABB{mi,ni3,q) 






ni2 


-'HABcimi,q, ni3) 






'HAAB{mi,m2,q) 






ni3 


-'HAAc{mi,m2,q) 


C 


mi 


Ti- ABC {m2, 1713, q) 




C 


nil 


'HACc{m2,q,i7i3) 




ni2 


-'HABc{mi,ni3,q) 






ni2 


'HACcimi,ni3,q) 




ni3 


-'HAAc{mi,ni2,q) 






ms 


'HAAc{mi,m2,q) 
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ABB 


'HABB{mi,m2,m3) 




ABd 


'HABdimi,m2) 


A 

/d 


mi 

7712 

ms 


riABBiq, m2,m3) 
Ti-AABiq, mi, ma) 
'HAAB{mi,q, m2) 




A 


mi 
m2 
d 


nABd{q,m2) 


'HAAB{mi,q,m2) 


R 


/ / fc]^ 

7712 
7713 


u 

'HABB{mi,q, ma) 
'HABB{mi,m2,q) 




R 
u 


Til -\ 

m2 
d 


T~(-BBd{m2, q) 

'HABd{mi,q) 

'HABB{'mi,m2,q) 


c 


mi 
m2 
w-3 


'HBBc{'m2,m3,q) 
'HABc{'mi,m3,q) 
-'HABcimi,m2,q) 




c 


mi 
m2 
d 


-'HBCd{m2,q) 
-'HACd{mi,q) 
-'HABcimi,m2,q) 














ACd 


T-(-ACd{mL, m) 




BBd 




A 


mi 
m2 
d 


'HACd{q,m2) 


'HACd{mi,q) 




A 


mi 
m2 
d 


-n-ABdiq, ^2) 
-'HABd{ini,q) 
-'HABB{q,rni,m2) 


B 


mi 
m2 
d 


-'HBCd{q,m2) 
-'HABdimi,q) 
-'HABc{rni,q,m2) 




B 


mi 
m2 
d 


' ^BBdVi-) '"'2} 

'H.BBd{mi,q) 



C 


mi 
m2 
d 




'HACd{mi,q) 
'HACc(,rni,m2,q) 




C 


mi 
m2 
d 


-'HBCdi'm2,q) 
-'HBCdimi,q) 
-'HBBc{mi,m2,q) 



ABC 


'HABc{mi,m2,m3) 


A 


mi 


T-i-ABciq, 'm2,m3) 




m2 


'HAAc{q,mi,m3) 




ms 


'HAAB{q,mi,m2) 


B 


mi 


'HBBciq,m2,m3) 




m2 


'HABc{'mi,q,m3) 




m3 


'HABB{mi,q,m2) 


C 


mi 


'HBCc{'m2,q,m3) 




m2 


'HACc{mi,q,m3) 




m3 


'HABc{mi,m2,q) 


d 


mi 


-'HBCd{m2,m3) 




m2 


-'HACd{mi,m3) 




ms 


-'HABd{fni,m2) 



6.2.4 Beyond Two Fermions 

Using the same methods as above and assisted by SWI's, the analysis can be 
extended to other external particles than gluons and two fermions. For instance, 

the ratio of coefficients between a {g~ , g~ , Sab, Sah) and a {g~ , g~ , g~ , g^) ampli- 
tude is found by squaring the ratio between coefficients of {g~ ,g~ , f~ , f~^) and a 



82 



,9 ,9 ,9^), e.g. 

c^"'i{m^,m^},{ml-},{q'+},d+) 



'HAAB{mi,m2,m3 
( -'HAAc{mi,m2,q)\ ^ 3^ 



K'HAABi'mi, 



c^-({mr, m2-}, {m^}, {q+}, d+). (6.35) 



If we moved four units of —1/2 helicity from 777,3 to q we would have to take 
the parenthesis to the power 4, which is of course equivalent to exchanging 
WaabI^t-Ii "72, 777,3)^ with WaacI^t-I) '"T'2, 9)^ in (|6.2ip as it should be. 

Amplitudes with both scalars and fermions can be constructed in similar 
ways by viewing them as pure-glue amplitudes where units of —1/2 helicity have 
been moved around. For each move, one of the original 7ik factors is removed 
and replaced by another which captures where the helicity was moved to. In 
general this does not fix the overall sign unambiguously, so we must resort to 
SWIs to be completely sure. As an example we can try to compute the coeffi- 
cient c'^™ ({ 777 7772 ^^^}' {"^3 ^^^}) {^'^^li '^^) t>y noting that it comes from moving 
— 1/2 helicity from 7772 to q in c^™'({mj~, Tra^}, {7/73 "^^^}, {g^/^}, d'*'). This tells us 
that the relative factor must be ±Ti.ABc{^i,n^3, Q)/'HAABimi,m2,m3), and that 
we must use a SWI of the type (j6.25p which describes moving helicity from some- 
where to in an AAB configuration to corner C . Such an SWI is represented in 
the first table above in the last three lines, namely 

= {mir])TCABc{'m2,m3, q) - (777-2??) Wabc ("J-i, 

-{m3r])7iAAc{mi,m2, q) - {qr])7iAAB{rni,m2,m3). (6.36) 

This leads us to conclude that 

c3-({7nr, m-'/'}, {7773^^/^}, d+) 

= ;''T^"^^'"^^'^i c3-({777r,777,-},{7773-^/^},{gV^},d+) 

T-CAAB{mi,'m2,m3) 



Ti-ABC ("T-i , m3 , q)TiAAC (m-i , 7772 , 1 



-c 



,3m 



({"T-l ,^2 },{m3 },{g+},d+). 



Ti- AAB (mi, 7712, nisy 

(6.37) 

With these methods, three-mass box function coefficients for any NMHV external 
particle configuration can be calculated. 



6.2.5 Beyond Three-Mass Boxes 

Up until now we completely focussed on three-mass boxes because they were 
MHV constructible and have a well controlled structure. Two-mass-hard boxes 
are also MHV constructible but they have two contributions for each coefficient 
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since either of the two massless corners can be the googly MHV one. We can eas- 
ily obtain those coefficients, however, from the three-mass coefficients by taking 
massive corners to have only one momentum on them, 

c^'^^iA, B, c, d) = c^'^iA, B, {c}, d) + c3"*({d}. A, B, c). (6.38) 

The two-mass-easy boxes are not MHV constructible because the opposite mass- 
less corners are both googly MHV and the two massive corners must be MHV and 
NMHV, respectively. Rather than using explicit expressions for the the NMHV 
amplitudes, we can use the same procedure as [37] where the known structure of 
IR divergences in A/" = 4 super- Yang-Mills tells us that 



^1— loop 



= -5E(^)■^"^^ (6-39) 



which again translates to linear conditions on the box coefficients. These can be 
inverted (at least for an odd number of external particles) to produce expressions 
for the two-mass-easy coefficients in terms of the three-mass coefficeints. Because 
we are considering = 4 supersymmetry in both external and internal states, 
these identities hold irrespective of the external states. Formally, the identity is 

c2--(A, b,C,d) = J2c'"'iC,...,...,b) + J2c'"'iA,...,...,d), (6.40) 

dec beA 

where the summations rum over all three-mass box coefficients that fulfill the 
condition. 

Just as two-mass-hard boxes were degenerate three-mass boxes, one-mass 
boxes can be characterized as degenerate two-mass-easy and three-mass boxes. 
This gives us the formula for the one-mass coefficients, 

c^'"(A, b, c, d) = c^'^^A, b, {c}, d) + c^'^iid}, A, {b}, c). (6.41) 

Since there are no further contributions to the one-loop amplitude, this consti- 
tutes a complete characterization of the one-loop NMHV amphtude in A/" = 4 
super- Yang-Mills. 
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Chapter 7 

One-Loop Amplitudes in AT = 8 
Supergravity 



This chapter introduces the calculation of one-loop amplitudes in M = 8 super- 
gravity and the "No Triangle Hypothesis''^ for the structure. It is based on the 
article [50] , which extends the results of earlier work [48] . 

7.1 One-Loop Structure of Maximal Supergravity 

The results of sections [5 . 1 . 1 1 and 1 5 . 1 . 2 1 apply just as well to a gravity theory as to 
Yang-Mills. Thus, any one-loop amplitude in four dimensions can be written as 
a sum of boxes, triangles, bubbles and a rational function. 



There are, however, no theorems that use supersymmetry to constrain this further 
as with Yang-Mills, and generically, any theory of gravity should be expected to 
include all types of terms. 

For a low number of external gravitons, there are some results that do con- 
strain the form. If we go back to the integrand of the loop integral, it will generally 
have m < n propagator terms where n is the number of external particles. Since 
all Feynman vertices in gravity contain momentum squared, we should have 2m 
powers of loop momentum in the numerator. By multiplying and dividing by 
all propagators not in the expressions, this can be written as 2n powers of loop 
momenum in the numerator and n propagators. 



^ Given the present weight of evidence it should be termed the "No Triangle Conjecture", but 
the original name seems to have stuck with it. Considering what it actually says, it ought to be 
called the "Box Hypothesis" or something similar. 






(7.2) 
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When we do Passarino-Veltman reduction on this, each propagator can be ex- 
changed for one or two fewer powers of loop momentum. For any n, this is 
expected to result in a combination of all the terms appearing in ()7.1|) . How- 
ever, we should also expect that there are supersymmetry cancellations like in 
Yang-Mills that reduce the power 2n. This is indeed the case in the so-called 
string-based approach to loop calculations [43, 80] (for the string-based approach 
in Yang-Mills, consult [27,42,45,46]) where cancellations can be made explicit 
which reduce the numerator power of loop momentum by 8 in = 8, 



This explicitly tells us that four-point amplitudes (as explicitly shown by a classic 
string calculation [94]) can only have a box integral, five-point can have boxes 
and triangles, six-point can have boxes, triangles and bubbles, while seven-point 
and above is again unconstrained. 

Somewhat surprisingly, the situation seems to be even better. It has been 
known since the end of the 90'es [24] that five-point and six-point MHV 1-loop 
amplitudes consisted only of boxes and it was conjectured that this would hold 
for all MHV amplitudes. When it was recently found that the box function 
part of the NMHV six-graviton amplitude was sufficient to explain all IR, soft 
and collinear properties of the amplitude, it was hypothesized that all one-loop 
amplitudes in AA = 8 supergravity consisted of boxes only. In particular, since 
there this means no triangles, it was called the No- Triangle Hypothesis. 

An equivalent way to express this at one-loop is that M = 8 supergravity 
has the same structure as AA = 4 super Yang-Mills, in particular that there 
are cancellations which remove additional powers of loop momentum from the 
numerator of integrands to reduce it to 



This is n — 4 more cancellations than M = 8 supersymmetry is known to give and 
the origin is unknown. As described in the introduction, this and many other 
results have led to speculations that M = 8 supergravity is perturbatively UV 
finite to all orders and that this may be caused by an unknown symmetry present 
in all gravitational theories, independent of supersymmetry, which softens the UV 
behaviour. 

7.2 Evidence for the No-Triangle Hypothesis 

The present strategy for supporting the No- Triangle Hypothesis is to determine 
the part of all six and seven-point amplitudes depending on box integrals and 
to argue that no further cut-constructible parts can be present. This leaves 
rational parts which are known to be absent at six-points; additional arguments 
are presented why these ought to vanish at seven-point and higher. 




(7.3) 




(7.4) 
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7.2.1 Box CoefRcients and Soft Divergences 

Box coefficients can be determined by use of the quadruple cut method of section 
15.41 There are two strategies for determining these; the direct and the KLT- 
inspired. The latter works by using the KLT relations on all corner amplitudes 
of the quadruple cut to write the M = 8 coefficient as a product of kinematic 
invariants and two = 4 coefficients [28]. As an example, the coefficient of 
/4(1~,2~,3+, (4+,5+)) given by the quadruple cut 



2" 3+ 




which can be KLT'ed to 

cf-t-,3+,(4+,5+) = ^A^(-/^l-,^2+)-M(-/2",2-,/+) 

xAl(-/3-, 3+, 4+, 5+, ) 

= —A.{—li,l ,12) •^{~^2 '^3^) -^(~^3 5 ) 

XS45^(-ZJ, 4+, 5+, l^)A{-q,5+,4+, 1^) 

= 2s45C^^|-,3+,(4+,5+)'^,2-,3+,(5+,4+)- C''-^) 

In the cases where there is a summation over the multiplet running in the loop, 
the rewriting will still hold because the M = 8 multiplet is two copies of the 
M = A multiplet, each being summed over in the two M = A coefficients on the 
right hand side. Unfortunately, this strategy becomes increasingly cumbersome 
at higher points. The KLT relations become longer and involve permutations of 
legs which lead to the N = A coefficients being the non-planar ones. 

The more direct strategy is to use the expressions for = 8 supergravity 
amplitudes. These can be derived from on-shell recursion in more compact form 
than those given by the KLT relations, and presumably give rise to the most 
compact forms of the coefficients. The coefficients for general MHV amplitudes 
were given in [24] and the coefficients for the six-point NMHV amplitude were 
given in [28]. Thus, the remaining seven-point amplitudes are the NMHV ones. 
A large portion of these are MHV constructible (c/. section EH]) but some one- 
mass coefficients may depend on the NMHV six-graviton tree amplitude. One 
even requires summation over the full multiplet, and we thus need to extend 
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the NMHV six-point tree amplitude to the case where two external particles are 
non-gravitons. This can be achieved by redoing the on-shell recursive calculation 
of [61] and the result is given in appendix C of [50]. The procedure like the one 
used in section 16.2.51 cannot be used as it relies on the previous knowledge that 
the amplitude contains only boxes, the fact that we are trying to establish. All 
the seven-point NMHV box coefficients are collected in appendix A of [50]. 

The results for the box coefficients constrain the rest of the amplitude because 
it must have the soft (1/e) divergence [81], 



_y^l-loop 



(4vr) 



Ei<jSijlog[-Si 



2e 



M'""". (7.6) 



If this holds for just the box integral part of the amplitude, the triangles and 
bubbles must have their soft divergences add up to zero. That this is indeed 
the case has been checked numerically for the all amplitudes up to and including 
seven-point. In practise, this is done by extracting the soft parts of all the box 
functions and rephrasing (|7.6|) as identities among the box coefficients. 

The soft divergences of one and two-mass triangles are best studied by rewrit- 
ing them in a basis consisting of the functions 

where runs over all independent kinematic invariants. Since the soft parts 
of this basis are independent, the vanishing of their sum amounts to requiring 
the vanishing of every coefficient. Going back to the basis of one and two-mass 
triangles this means that all coefficents of these are also zero. The three-mass 
triangles have no IR singularities and can thus not be ruled out in this way. The 
bubble functions are proportional to 

= I - log(-i^^) + . . . , (7.8) 
so all we can say about them is that their coefficients add up to zero. 



7.2.2 Ruling Out Bubbles 

The safest way of ruling out bubble contributions is of course to calculate them 
directly. For this, the article uses two methods which are closely related. The first 
is the method of Britto, Buchbinder, Cachazo and Feng [58] where the unitarity 
cuts of amplitudes are explicitly integrated and the results compared with the 
cut of the expansion in terms of scalar integrals. By inspecting all unitarity cuts 
of seven-point NMHV amplitudes in this light, bubble integrals can indeed be 
ruled out. 

Rather than going through this technical, albeit interesting, calculation, we 
will use a method which ties together the behaviour of tree amplitudes deformed 
by a complex parameter z as in on-shell recursion (c/. section [3^ and the UV 



88 



behaviour of one-loop amplitude. This choice is made to illustrate one of the cen- 
tral points of this thesis, namely that properties of graviton tree amplitudes (such 
as z — > oo behaviour) are tightly linked with the whole perturbative expansion 
of = 8 supergravity. 

Using the expansion ()7.ip we can write the unitarity cut in a particular channel 
with momentum P and lift the momentum integration. This gives 

where C and V are the sets of box and triangle functions which have a cut in the 
relevant channel and the / — -fCi.2,3 represent the remaining propagators. Only 
one bubble function has the cut; this is represented by e^. Finally, we have added 
a total derivative term which is set to zero by the integration. By the normal 
logic of unitarity, this integrand should be equal to a product of two on-shell tree 
amplitudes, summed over internal helicities, 

^ (7.10) 

/ll,/l2=± 

We can now deform this product in the same way as done in the BCFW con- 
struction of on-shell recursion described in section 13. 3t namely by making the 
shift 

\h) = \h)+z\h), \h] = \h]-z\h]. (7.11) 

Under this shift, the integrand of boxes and triangles go as z"^ or as z ^ co 
and bubbles go as z". This immediately suggests a condition for absence of a 
bubble: If (|7.1U|) tends to zero as z ^ oo (and D{li, I2) does as well) the bubble 
corresponding to the cut is absent. This condition can be related directly to an 
equivalent condition for vanishing of bubbles in the direct integration method 
mentioned above, so it is indeed an appropriate check. 

The behaviour of tree amplitudes under the under above shift ()7.1ip is not 
known in general, but two things are known: For MHV amplitudes where the 
helicities of (/i,/2) are given as (/ii,/i2)i then the scaling can be derived from 
known expressions for the amplitude [12, 17] to be 

^^^^ (/,^,/,2) = ( + ,+),(_,_),(+,_), 
_y^MHV_^6 ^YlQTi (/li,/l2) = (-,+), 

^MHV _ ^2h+2 ^Yien {hi, h2) = {-h, h). (7.12) 

It is also known generally [14] that 

Mr^z-^ when (/i^, /la) = (+, -) (7.13) 

for all graviton amplitudes. Furthermore it can be verified that the two first lines 
of ()7.12p holds for all six and seven-point tree amplitudes needed here. Thus, it 
is not unreasonable to assume that this is a general picture. 
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The known z — > oo behaviour of MHV amphtudes can be used immediately 
to derive the behaviour of the cut integrand for one- loop MHV amplitudes. A 
singlet cut integrand where there is only one state running in the loop has the 
generic form 

M{l-,2-,3+, . . .,r+,lt, -l+)M{-q,q,{r + 1)+, . . . ,n+). (7.14) 




Under the shift, both tree amplitudes go as so the cut goes as and there 
are no bubble contributions. 

If the cut is non-singlet, the whole multiplet can run in the loop, 




the integrand becomes 

' 14-2.., 

xA^ (2-, 3+, . . . , r+, l^^^ (r + 1)+, . . . , n+) (7.15) 

However, we know from SWIs that MHV amplitudes with different external par- 
ticle content are related, e.g., 

A^(2-,3+,...,r+,/^-/2')= (my ''A^(2-,3+,...,r+,/+ -/2-), 

(7.16) 

so we can write (I7.15P as 

p X 7W(2',3+, . . . ,r+,Z+, -q)M{-l^,lt, (r + 1)+, . . . ,n+) (7.17) 
where 
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Under the shift, the two tree ampUtudes go as and z^, respectively, while p, 
whose form is a consequence of AA = 8 supersymmetry, goes as z^^. In total, the 
integrand goes as z~'^ which rules out bubbles in all one-loop MHV amplitudes. 

For singlet NMHV cuts, the situation is quite the same as for MHV; the 
integrand is 

M{lY,...,-l^)M{lt,...,-lt), (7.19) 

and since both tree amplitudes go as z~'^ under the shift, the integrand goes as 
z~^. This explicitly holds true for all cuts of six and seven-graviton amplitudes, 
and continues to hold for any singlet cut if the two first lines of (|7.12|) hold beyond 
MHV. 

The non-singlet NMHV cuts are more involved since they require the sum- 
mation over the M = 8 multiplet. For some integrands, such as 

E(-l)''(4_^/,)-^(l"'2"'4+,/r\-^2)-MH?,^2"',3-,5+,6+),(7.20) 
h ^ ^ 

The calculation is essentially the same as for the MHV case because the left 
NMHV amplitude is also a googly MHV amplitude where 

Al(l-,2-,4+,/r',-/2')= (-|^)' ''A^(l-,2-,4+,/r,-/2+), (7.21) 
such that the integrand becomes 

( (Sm ) '-^(1"' -^2+)-^(-^^ 3-, 5+, 6+), (7.22) 

and goes as z~'^. But in the seven-point calculation we will need cuts like 




with the integrand 

E(-l)''(4_^J-^(l"'2"'4+,5+,/r\-^2')-M(-^^/2"\3-,6+,7+), 
h ^ ^ 

(7.23) 

the expression depends on the six-point NMHV tree amplitude. As mentioned 
above, this has been calculated in the paper and is found to contain 14 terms. 
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of which 12 have the structure X x y2/i-4 ^-^q last two only contribute at 
/i = ±2. When this expression is inserted into (|7.23p we find that it goes to 
zero as z ^ oo term by term (after the h summation) although it does not go 
as well as z~^. Again this highlights the connection between recursion relations 
and UV structure: Using recursion relations we can write down a form of the six- 
point amplitude which rules out certain bubble contributions to the seven-point 
one- loop amplitude without the need for further cancellations. 

Since we have now calculated the z — > oo behaviour of every independent 
cut at six and seven-points, and found that they approach zero, we conclude 
that there are no bubble contributions to the one-loop amplitudes, and that the 
absence can be tied to recursion relations for gravity. 

7.2.3 Ruling Out Three-Mass Triangles 

As mentioned above, three- mass triangles have no IR singularities and have to 
be ruled out explicitly. This can be done by inspecting the triple cuts 

hi,h2,hi 

xM{..., -l^''\l^'^)M{. . . , -q''\l'^^)M{. . . , (7.24) 
and arguing that they are all results of triple cuts of the boxes, 

/ d^'hSilDSilDSill) n (7-25) 

where li — K{i) is the momentum of the remaining uncut propagator. The helpful 
fact here is that there are kinematic regimes where the three-mass triple cut is 
real, and thus the remaining integration can be performed numerically without 
appealing to complex Minkowski space methods. 

The box functions that have three-mass triple cuts are the 2-mass hard, 3- 
mass and 4-mass boxes. Since we are only going up to seven-point amplitudes 
we need not worry about the last of these. The result of the cut integrations are 



h(ki,k2,Ks,K4,) 
h{KuK2,K3) 



vr 



2{h + k2)^{k2 + Ks)' 
IT 



2 (fci + ^2)2(^2 + K3)2 - KiKi 



vr 



cut 



2y'Kt + Kl + Kl - 2{KiKl + Km + KlKl) 

(7.26) 

With the knowledge of the box coefficients we have, it is fairly straightforward 
to numerically integrate all the triple cuts and compare them to the expectation 
from the boxes at some randomly chosen kinematic points. This explicitly rules 



92 



out three-mass triangles for both six and seven-point amphtudes. For the six- 
point NMHV, the independent triangle coefficients are 



d3-({l-,2-},{3-,4+},{5+,6+}) 

d3-({l-,4+},{2-,5+},{3-,6+}). (7.27) 

The triple cut corresponding to the latter is a non-singlet cut which requires a 
sum over the M = 8 multiplet. This means that it is non-zero for M < 8. The first 
triangle coefficient can be expected to be zero in any massless theory of gravity. 
For the seven-point NMHV, the independent three-mass triangle coefficients are 

d3-({l-,2-},{3-,4+},{5+,6+,7+}) 
d3-({l-,2-},{3-,4+,5+},{6+,7+}) 
d3™({l-,2-,4+},{3-,5+},{6+,7+}) 

^^^({l-, 4+}, {2-, 5+}, {3-, 6+, 7+}). (7.28) 

Again, the three first correspond to singlet cuts which are only sensitive to the 
presence of the graviton and should be zero in any massless gravity theory, while 
the last is non-zero for M < 8. 



7.2.4 Factorization and Rational Terms 

We have now ruled out contributions to the cut-constructible parts from anything 
but boxes at six and seven-point. At six-point it is already known that there are 
no rational parts, so this provides a verification of the No- Triangle Hypothesis 
in that case. For seven-point amplitudes, however, there is no general argument 
ruling out rational contributions. 

Computing rational contributions is in general very hard in gravity theories. 
A strategy that works in Yang-Mills theory is to do the full calculation in 4 — 2e 
dimensions (c/. section I5.2.3P and extract the rational pieces by taking e ^ 0. 
This works well because any massless field content can be decomposed into an 
AA = 4 multiplet, some number of AT = 1 chiral multiplets and some number 
of scalars, and because it can be proven explicitly that the two supersymmetric 
theories have no rational parts; only the scalar loop calculation needs to be done 
exactly in 4 — 2e dimensions. For gravity, there is no such theorem saying that 
supersymmetric theories have no rational contributions, so the calculation needs 
to be done in 4 — 2e dimensions for all states of the multiplet, thereby rendering 
most of the benefits of the spinor-helicity formalism useless. 

An alternative is to use the methods of section 15.51 or, more generally, the 
known factorization structure of one- loop gravity amplitudes. As in the gauge 
case, this means in particular multiparticle, collinear and soft factorization. The 
multiparticle factorization is, on the grounds of general field theory arguments 
[30] , the same as for gauge theory (15. 7p with the obvious corrections for the lack of 
colour ordering in gravity. The collinear behaviour at one-loop is actually simpler 
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than in Yang-Mills (15. 8p as the splitting functions have no loop corrections. This 
gives [24] 

>^l-'°°P(a'^^6'^^...) ^ ^ Splitl7(z,a'^»,6^'')7Wi-i°°P(P'^,...) (7.29) 

h=±2 

as a ^ zP, b ^ {1 — z)P. The splitting functions for gravitons are 
Splitf^""(z,a+,6+) = 



Splitr"(^,a+,6^ 



1 [ab] 
z{l — z) {ab) 



Splitr^(z,a-,6+) = -T^j^. (7.30) 

1 — z [ab) 

There is also a universal soft behaviour when a graviton momentum approaches 
zero, which also does not incur loop corrections, 

A^^'°°P(.±, . . .) - Sr^{s^)Ml-yr{- ■ •), (7.31) 

as kg 0. When the graviton going soft is number n in some arbitrary ordering, 
the soft factor is 

^~ (ln)(n(n-l)) ^ (in) ' ^''^ ' 

These limits put some quite strong constraints on the possible rational parts 
which must have the same factorization as the cut-constructible parts. We know 
that the rational parts of the six-graviton amplitude are zero, so by ()7.29p and 
(j7.3ip the rational part of the seven-graviton amplitude can have no collinear or 
soft singularities. It is hard to imagine any rational spinor expression which has 
none of these singularities while being of mass dimension 4 and having the right 
holomorphic weight of spinor s as determined by their helicities. Indeed, it is hard 
to imagine such an expression for any number of external gravitons, so by "loose 
induction" we should expect that all rational parts are zero. 

This argument could probably be formalized by using on-shell recursive meth- 
ods like those of section 15.51 for gravity. Some progress has been made in this 
direction [53] but the approach is still too undeveloped to be used for the problem 
at hand. 



7.3 Conclusions 

Using an array of different methods, we have explicitly calculated the cut-containing 
parts of the NMHV six and seven graviton one-loop amplitudes in M = 8 super- 
gravity and found that only box integrals contribute. At six points, this deter- 
mines the amplitude completely while at seven-point we have argued that an 
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additional rational term is highly unlikely. This proves the No- Triangle Hypoth- 
esis at six points and renders it highly credible at seven. 

The methods used can be applied at any number of external particles, and it 
does not appear that they would suddenly start failing. This applies in particular 
to the ruling out of bubbles which relied on very general scaling arguments for the 
amplitudes, scalings that arc expected to hold to any number of points. Together 
with the arguments that the rational parts of amplitudes are likely to be zero, 
this provides strong evidence that the one-loop amplitude in = 8 supergravity 
consists exclusively of boxes just as A/" = 4 SYM. 

As explained in the beginning of the chapter, this amounts to the cancellation 
of n — 4 powers of loop momentum in a n-point amplitude, a cancellation which is 
unexplained. Similar cancellations seem to occur at higher loop order also. They 
have been explicitly found in the four-point amplitude at two [23] and three [25] 
loops, where it was proven that the divergences of maximal supergravity and 
maximal super- Yang-Mills are the same in any dimension. This gives a strong 
hint that the divergences of the two theories match in general, and that M = 8 
supergravity in four dimensions is UV finite in particlular. It also hints that 
there is some unknown symmetry or dynamical principle present in A/" = 8 which 
is responsible, although there are no candidates for such a symmetry. There is 
evidence that the cancellations occur already in non-super symmetric gravity [21], 
and thus the symmetry should be visible in some form without the presence 
of supersymmetry. The identification of such a symmetry will be of immense 
interest. 
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Chapter 8 

Amplitude with a Higgs 



In this chapter we make the full determination of a particular one-loop helicity 
amplitude, ^^(Higgs, 1^, 2~, 3^^, 4+). The chapter is based on [10] but presents 
the calculation in a somewhat different way. The article calculates the cut con- 
structible parts, using the one-loop MHV rules of section I5.3|, of the all-n MHV 
amplitude ^^((/), 1~, 2~, S"*", . . . , n+) (0 to be defined) and the rational part of 
„4^(Higgs, 1~, 2~, 3"*^, 4''") with on-shell recursion from section [^31 In this chap- 
ter, we limit ourselves to the four-gluon case all the way through and compute 
the cut-constructible parts using classical unitarity as in section I5.2i 

Even though the pp — > Hjj process, where this amplitude is relevant, is 
not included in table 11.11 of wanted NLO calculations, it is generally considered 
of the same importance. Some semi- numerical results already exist [84], but 
clearer analytic results are still missing. The process forms a background to 
Higgs production by vector boson fusion [65]. 

8.1 Higgs in the Large Top Mass Limit 

In the present calculation we will make the assumption that the top mass is large 
or, to be exact, that the kinematic scales involved in the scattering are small 
compared to twice the top mass. The dominant coupling of the Higgs to massless 
QCD particles comes through a top quark triangle loop connecting to two gluons, 



In the large top mass limit, such an interaction is well approximated by an effec- 
tive term in the Lagrangian [130, 132] 




2 ^ 
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where H is the Higgs field and C is an effective coupling constant equal to [97] 
as A , lla^ , 2^ 



6nv \ An 

In the addition to the Higgs, we can add a pseudoscalar Higgs A such that 
the total effective Lagrangian becomes 

^ [hTiF^^F^"' + z^TrF^/F'^^) (8.3) 

where 

If we now change variables to 

(/)=i(// + M), </>t = l(i7_M), 

FsD = liF"" + *F'^''), F^j, = ^{Ff^-^ - ^Ff^n, (8.5) 
the effective Lagrangian becomes 

C [<l,l^FsD^.uF^sD + <P^'^^FASD^..F^^'sD) > (8-6) 

and Higgs amplitudes can be recovered by 

A{H,...)=A{<^,...)+A{(t)'^,...). (8.7) 

The Higgs is colourless and does not participate in the colour ordering, but is 
conventionally written first. 

The central observation is that tree amplitudes for 4> and (j)^ are separately 
simpler than amplitudes for H ox A. In particular, they follow a pattern like that 
of the Parke-Taylor amplitudes ()2.4ip , 

A\^, l+,...,n+)=0 
^°((^,l~,2+,...,n+) =0 



^"(</.,l~,2+,...,r,...,n 
for 6, and 



(12)(23)---(nl) 



^O(0t,l-,...,n-) = O 
^0(,^t^l+,2-,...,n-) = 

^0(^t,i.,2-,...,.-,...,n-) = (-ir^J^, (8.9) 

for (f)^ . Even though the amplitudes look like the MHV and googly-MHV, the con- 
dition between the gluon momenta for conservation of momentum is different be- 
cause of the presence of the cj) momentum. Notice also that parity changes the sign 
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of gluon helicities and daggers the (p, so there are no easy +, +, — , .) 
ampHtudes. In fact, there is a simple (;/)-an-minus amphtude (and 0^-all-plus 
ampHtude) , 



4 



[12] [23] •••[«!]' 

4 

^°<^^-'"--"")= (12X23) ..(nl) - e*'!") 

The similarity of the (j)-MHV amplitude to the MHV amplitude have led 
people to develop MHV rules for amplitudes with a single 4>, where a 0-MHV 
vertex must be used once but where the off-shell prescription is unchanged [8,77]. 
The extension to one-loop amplitudes was explored in [6] and the article on which 
this chapter is based [10]. The former deals with the cut-constructible part of 
A^cl^, -,-,...,-). 

One-loop amplitudes involving a (j) and at most one negative helicity gluon are 
purely rational because the tree amplitude vanishes. All-n expressions for these 
amplitudes have been constructed in [20] where Schmidt's results for a Higgs and 
three gluons [129] is converted into the 0/(/>^ notation. Some of those we need 
are given at the relevant point below. 



8.2 A One-Loop Higgs Amplitude 

The methods presented in chapter can be applied just as well to amplitudes 
containing a Higgs in the large mt limit. Here we will apply them to the cal- 
culation of A{H, 1^, 2^, 3^", 4+), that is, the one-loop amplitude for a Higgs and 
four gluons with the 'adjacent MHV helicity configuration. This configuration 
has a non-zero tree amplitude and we can thus expect cut-containing pieces at 
one-loop. In practice we will calculate A{4>, 1~, 2~, 3"'', 4^) and use that 

A{H, 1-, 2-, 3+, 4+) = Ai^, 1-, 2-, 3+, 4+) + ^(0^, 1", 2", 3+, 4+) 

= ^(0, 1-, 2-, 3+, 4+) + ^(0, 3-, 4-, 1+, 2+)*. (8.11) 

^(</), 1^, 2^, 3+, 4+) is split into its cut-containing pieces and its rational pieces. 
For generality we calculate the cut-containing pieces for A{(l), 1^, 2~, 3^, . . . , n^). 
The whole calculation is done in the four-dimensional helicity (FDH) renormal- 
ization scheme. 



8.2.1 Cut-Constructible Pieces 

As described in section 15.1.31 the states running in the loop can be written as 
a linear combination of an = 4 multiplet, an M = 1 chiral multiplet, and a 
complex scalar referred to as AA = 0. The calculation of cut-containing pieces can 
be performed separately for the three multiplets, and the result assembled from 
them. 



99 



Our plan of action is to first uncover any box functions by the use of quadru- 
ple cuts. We then consider the (two-particle) unitarity cuts where we isolate 
the terms coming from the (by now known) boxes and relate the remainder to 
triangles and bubbles. 



Quadruple Cuts 



If the amplitude is cut in four places, the four corners must together have six 
minuses (two external and four internal). These must necessarily be distributed 
as two three-point googly corners (opposite) and two MHV corners. Since the 
corner with the cannot be three-point googly, it must be one of the MHV 
corners. For = 4, which contains the gluons, this gives us the following six 
types of non-zero quadruple cuts: 





Notice that all permitted cuts necessarily have gluons running in the loop. This 
immediately shows us that M = 1 and = cannot contribute since they do 
not have gluons running in loops. 

The calculation of the prefactors of the box functions is almost identical to the 
one without the (j). Thus, we will only present one example calculation, namely 
that of the coefficient c of the box corresponding to the second cut above. From 
(|B.3|1 in appendix [B] we derive that 

,m— 1 1 

t1 — ; -T , t2 



(m2) (m2) 

"T-f3,m~l| , _ \m){21 

(m2) ' ^ ~ (m2) 



^ ^ |2)(mP3,^^i| ^ ^ |m)(2P3,^_i| 



such that 

^ 1 OM! M-h)f 

2 {l2{-h)){{-h){m + l))((m + 1) • • • 1) [(-/2)2][2/3] 

^ U-h)? [{-k)mf 

((-«3)3)(3 • • • (m - l))((m - l)U) [mh][h{-h)] 
^ 1 (12)(23)((m - l)m)(m(m + l))(l/2)^[^2/3]^(^3^4)^[/4m]=^ 

" 2 (• • •)(/2/l)(/l(m + l))[/22][2/3](/33)((m - l)/4)Mi][/i/4] 

^ 1 (12)(23)((m - l)m){m{m + l))(12)3(-(mP2,m-iP3,m-im))3 
2 (• • •)(m2)4(2m)(m(m + l))(mP2,m-i2](-[2P3,m-im))(23) 

^ (2m)3(2P3,^-im]3 

((m - l)m)(-[mP3,„2))(-(2P3,^P3,^_i2)) 
(12) Y (2P3,™_im](mP3,„_i2]\ 

-^tI ^ J- ^'-''^ 

As usual for MHV amplitudes, this is just the Gram determinant times the tree 
amplitude. This holds for all coefficients, so the boxy part of the amplitude is 
given by 

,^2\4 mi-l 

i+l,m2-l,'^2), (8-14) 

mi=lm2=mi+2 
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which is a combination of 1-mass and 2-mass-easy box functions. The difference 
in comparison to the result without the (p is primarily that there are twice as 
many 2-mass boxes because the cj) makes the two massive corners non-identical. 

Two Particle Cuts, = 4 Case 

By counting minuses, we can again see that the two tree amplitudes separated 
by the cut must be MHV and (/)-MHV, respectively: 




Regardless of the position of 1 and 2, the cut integrand becomes 

(12)^ ((mi - I)mi){hl2){m2{m2 + l)){l2h) 
(. . .) ((^^ _ I)h){hmi){m2l2){h{ni2 + 1)) ' 



.15) 



For generic values of mi and m2 it can be shown, analogously to the pure gluon 
case, that this is equivalent to the same cut of 



(12^ 
(...) 

+F4(mi,Pm,+i_m2-l,'^2) + F4 (mi , P^i +I,m2 , "l2 + 1)). (8.16) 



(^Fi{mi - l,Pmi,m2-i,"i2) +Fi{mi - l,Pmi,m2,rn2 + 1) 



The pure glue requirement that mi — 1 7^ m2 + 1 is, however, not present when 
there is a on the left tree amplitude where the cut 




must be taken into account. The decomposition is the same as above, but the 
box function F4(mi — 1 , ^^.2 ; "^2 + 1) = F^i^, Pm+i,m-iin^) is physically 
nonsensical as it stands. Luckily, its cut can be rewritten, 

1 ((mi - l)Pmi,m2("T'2 + 1) (mi - 1)] 

2 (Zi_(mi-l))2(;2 + (m2 + l))2 
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_ 1 {mPm+i^m-iinf 
2 {mlirn\{rnl2rn] 
_ {mPjn+i,m-im] {mPm+i,m-im] 

2(/i-m)2 + 2(/2 + m)2 ' ^ ' 

as a sum of the cuts of the two triangle functions. But since a triangle functions 
is independent of the ordering of its momenta, these two triangles are actually 
the same, leaving us exactly with the cut of ^3(771, Pm+i,m-i)- There are no 
contributions from bubble functions in the AA = 4 case. 

The conclusion is, that the amplitude should have the following dependence 
on triangles 



(12)^ 



Y^F3im,Pm+i,m-i). (8.18) 



Two Particle Cuts, J\f = I and = Cases 

With the two particle cuts, we have a possibility to let a full multiplet run in the 
loop if we have 1~ and 2~ sitting on different sides of the cut, that is, in the cases 




m + 1 




where m runs from 3 to n in both. For simplicity, we only treat the first of these. 
For AA = 1, the cut is given by 

(l/l)(l/2)(2/l)(2^2)( - ((l/i)(2/2) - (l/2)(2/i))') 
{hh){hini + l))((m + 1) • • • I){lh){l2h){h2){2 ■ ■ ■ m){ml2) 
^ (12)3(m(m + l)) {Il2){l22) 

(•••) (m/2)(/2(m + l)) 

_ {I2f f {ml){2l2) (l(m + l))(2/2) \ 
(mh) ^ ((m+l)/2) J 
_ (12)Y(2/2m](ml) (2/2(m + l)]((m + 1)1) \ 

(•••)V {mhm] ((m+l)/2(m+l)] ^' ^ 

We have now written this as cuts of vector triangles. To put this into the standard 
basis of scalar integrals we need to perform a Passarino-Veltman reduction [121]. 
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The details of this procedure is given in appendix \C\ and the result is (using 

P = —P2,m)- 

(12)Y(2P2,n^ml) (2P2.>n(m+l)l) X 

(• • •) V {'mP2 ] ((m + l)P2,m(m + !)];■ 

Here we have to remember that we have suppressed the integration, so what this 
really tells us is that the cuts we have calculated are the cuts of 

^ (12)3 ^ (2P2,.»ml) (2P2,n^(m + l)l) \^,^ 

{■■■)\{mP2,^m] ((m + l)P2,.„(m + l)]r^^^^'"^^- ^^'"'^ 



cr 

m=3 



Using that F2(-P2,m) - -?^2(-P2,m-i) = - log(P2^,m./-ff,m-i)' ^au rewrite this as 



- 2^ . o log ■ (8-22) 



The contribution from the second class of diagrams is 
(12)3/ (lP^_l(„^-l)2) {lPm,im2) 



cr 



m=3 
or 



(...) V((m-l)P^,i(m-l)] (mP™,H j^^iPm.ij, («.2^j 



(12)3 - (lP^,i(m-l)2) / ^^.1 \ 



The = case is complicated by the need for integral reduction. However, 
if we plunge right into it, the cut corresponding to the first diagram above is 

(l/l)2(l/2)2(2/i)2(2/2)2 



(/i^2)(^2(m + l))((m + 1) . . . l)(l/i)(/2/i)(/i2)(2 . . . m){ml2) 
(12)(m(m + l)) (l/i)(l/2)^(2/i)(2?2)' 
(...) {hl2f{ml2){h{m + l)) 
{12) {Ih){ll2){2h){2l2)'' ( {Im) (l(m + l)) 



(...) (/i/2)2 V("^^2) (/2(m+l)) 

, (12) (Hi/2l)(2/i/22) / (2/2(m + 1)1) _ (2/2ml) \ 
'(...) y^^rn+lMm+l)] {mhm] ) 

, (12) (lP2,>n^2l)(2P2,n^^22) / (2/2(m + 1)1) _ (2W)\ 

'(...) P|^ ^((^ + 1)^2(^+1)] {mhm]) 



(8.25) 



These are cuts of triangles with three-index tensors in the numerators. Again, 
we use the Passarino-Veltman reduction from appendix O (using P = —P2,m) to 
reduce 

(lP2,m?2l)(2P2,r«?22) (2?2ml) 

Pin. irnhm] ^'-"'^ 
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to 

(lP2,mml)(2P2.mm2)(2P2,mml) (lP2,m7^1)(2P2,>»m2)(27^ml) 
3(mP2,r„m]3 12(mP2,m"7,]2 
(1P2,^7^1) (2P2,^7m2) (2P2,mml) 

_ (lP2,^ml)(2P2,^m2)(2P2,^ml) ^ (2P2,^m2)(lP2,^ml)(12) 
3(mP2,mm]3 6(mP2,r„m]2 

(12)2(2P2,^ml) 
6(mP2,mm-] 

_ (lP2,^ml)(2P2,^m2)(2[P2,^,m]l) (12)^(2P2,^ml) 

6(mP2,mm]3 6(mP2,mH 

Thus, the bubble functions coming from the first diagram are 

(12) /(lP2,„^(m+ l)l)(2P2,„(m + l)2)(2[P2,„,(m+ 1)]1) 

Cr 



(8.27) 



3(---)V ((m + l)P2,„(m + l)]3 

(12)2(2P2,^(m+ 1)1) (lP2,^ml)(2P2,^m2)(2[P2,^,m]l) 
((m + l)P2,„^(m + 1)] (mP2,^m]3 

(mP2,^m] ) 
while those coming from the second are 

(12) / (lP„,i(m-l)l)(2P^,i(m-l)2)(2[P^,i,(m-l)]l) 



3(---)V ((m-l)P„,i(m-l)]3 

^ (12)^(lP^,i(m-l)2) ^ (lP^,iml)(2P^,im2)(2[P^,i,m]l) 
((m - l)Pm,i(m - 1)] (mP^,im]3 

Again, these can be rewritten as 



(12) ^ 

^ ' m=4 



/ (lP2,rnml) (2P2,,„m2) (2[P2,^, m] 1) 
V (mP2,mm]3 

(12)2(2P2,^ml)\ ( Plm 
'log ' 



+ 



(mP2,„,m] J V^2,m-i. 
(lP^,i(m - l)l)(2P^,i(m - l)2)(2[P^,i, (m - 1)]1) 
((m- l)P„,i(m- 1)]3 
{12f{lPm,iim-l)2)\^^^ ( Pl^i 



((m-l)P^,i(m-l)]y^"" V^'^-U 



(8.30) 
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Summing Up 

We can now sum up the cut containing pieces for arbitrary fermion and scalar 
content. In place of the logs we use 



Li{s,t) 
The result is 

(12)^ 
(...) 



log(g/t) 
s-t ' 



L3is,t) 



logis/t) 
(s-t) 



3 ■ 



^.31) 



{n , mi-1 s 

+l,nii-l) + F/^{m\^Pmx+\,m2-\i'^'l)\ 
ra\=\ ^ m2=mi+2 ^ 



N N 



+ E 

m=4 

(lP2,mml)(2P2,mm2)(2[P2,m,m]l) ^ . p2 ^ 

(lP^,i(m - l)l)(2P^,i(m - l)2)(2[P„,i, (m - 1)]1) 



3(12)3 



m— 1,1 . 



11 2nf Us 
' Y ~ 3iV ~ 3iV 



(2P2,mml) 
(12) 



(lP^,i(m-l)2) ^ 2 p2 



.32) 



We can now specialize to n = 4, since this is the amplitude we ultimately want 
to calculate. There, the result is 



^ ' m=l ^ 



(m, m + 1, m + 2) + F|'"^(m, P, 



m+l,m+2 



, m + 3) 



+F-(™,P„«,„«))+cr||(l-^ + | 



.33) 



(1P23 42)(2[P23,4]1)L3(S234,S23) + (2P4131)(2[P41,3]1)L3(S41,S341: 

[341(12)2 /^^ 2nf ns\/, , ^ , A 

+"^^34r U " W - 3iV J + ^i(«4i, .341) j . 



8.2.2 Rational Pieces 

To compute the remaining parts of the amplitude we use the methods described 
in section 15.51 The total one-loop amplitude is 



^1 = C + p, 



^.34) 
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but C and R contain unphysical poles which cancel against each other, but pre- 
vent us from finding R directly from recursion. Instead, we add and subtract a 
term CR with the same unphysical singularities as C, 



= {C + CR) + (R-CR), 



i.35) 



such that only physical poles are left in each of the parentheses. Our first ob- 
servation about (j8.34p will be that the unphysical poles only reside in the L3 
functions, and by rewriting it as 



L3(s,t) 



1 log 

2 {s - 
1 1 
2t(s-t)2 

1 /I 1 

2 V7 s 

-(- - 

2[t^ s 



1 + (s - t)/t 



t)3 
1 



+ 



1 



log 



l + {t-s)/s 



1 



4t2(s 

1 
1 



1 



(t-s) 

1 



2s{t 



+ 0{{s-tn 



we see that CR must be constructed by replacing L3{s,t) by 

1 



+ Oi{s-tf) 



i.36) 



-(- - 

2\t^ s 



i.37) 



Thus, 



CR 



cr[34]/ Hf ... 
6(34) V N N 
(2P4i31)(2[P4i,3]l 



+ 



(3P4i3]^ 



(1P23 42)(2[P23,4]1) 
(4P234]2 

1 1 ^ 
— + — 

S341 ■S41 



— + — 

■S234 ■S23 



^.38) 



Knowing now that R — CR has only physical poles, we can attack it with 
recursion relations. For reasons which will become clear later, we choose the 
recursive shift 



1 



11) + z|2) 



|2] = |2]-z|l]. 



^.39) 



Rather than proving that this shift has the proper 2; — >■ 00 behaviour, we will note 
that it does indeed work without the ip [89], and that the tests to be explained 
in the next section show no signs of bad z ^ 00 behaviour. The results of doing 
recursion on Cii are the overlap terms O23) O234, O41, and O341, coming from 
the singularities in S23) eic. These are found by leaving the kinematic invariant, 
shifting its coefficient according to (|8.39p . and inserting the z which puts us 
exactly at the pole. Taking O23 as an example, we note that CR has a pole 
where 



= S23 = (23) [32] = (23) ([32] - z[31]) 



[32] 
[31]' 



^.40) 
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such that 



|1) 



IA233] 
[13] ' 



|2] = |3] 



[12] 
[13]' 



23 



|A23l][3| 

[31] 



{8A1) 



and find that 



O23 



cr 
' 6 



-^(1 



6 

- _^ 
6 

The other overlaps are 



n/ ^ n,\ [34](TPi2342)(2(34 - 4Pi23)l) 

N n) (34)(4P234]2S23 

nf ns\ [34]si23[34](42)(si23(24)[43] - (23)[34](4Pi2 3]) 

(34)(4P231]2[34]2S23 



N ^ N J 



,\ [34]si23(42)(si23(42) + (4Pi23](32)) 



n f n c 

'N^NJ (34)(4P231]2S23 



(8.42) 



O234 = — 
O41 = - 



cr 
6 

cr 



nf Us \ [34]si234(42) (si234(42) + (4P123P342)) 

(34)(4P231]2S234 



N N I 



Uf [34](12)^ 
6 V N Nj (34)s4i ' 



__cr( _nf_ [34](2Pi34]((2Pi34] + (23)[34]) 

O341 -—Qy--lv + iV J (34)[41]2s34i 



(8.43) 
(8.44) 

(8.45) 



The hattings used in these and later calculations are given in the table below: 



i 


|1) 


|2] 


Pi 


23 
234 
41 
341 


|^'l233J 

[13] 

P1234P342) 


"'J[13] 
P34P1234I] 


|A23lJ[3| 

, [31] , 

Pl234lJ(2P34 


(2P34I] 

' ^ (24) 
P34P12342) 


(2P34I] 

1^4124) 

(24) 
P1234P34I] 


(2P34I] 
|4)(2P412| 

(24) 

-P34l](2Pl234 


(2P34II 


(2P34II 


(2P34II 



We can now turn to the results of doing recursion on R. This does not give 
us the R since we are summing over a limited set of poles, but rather gives 
us the 'direct recursive' terms DR23, etc, composed of (a sum of) an unshifted 
propagator times a shifted tree amplitude times a shifted rational part of a one- 
loop amplitude. One complication is that the three-point one-loop vertices cannot 
really be interpreted as the corresponding amplitudes, rather they have to be 
deduced from the one-loop splitting functions. In the case at hand, however, we 
are considering a scalar running in the loop and wc have chosen our shift such that 
the two external gluons in a three-point vertex have opposite signs. Fortunately, 
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the rational parts of the corresponding splitting functions are all zero (c/. ()5.39p ) 
so there are no contributions in our case from three-point one-loop vertices. The 
non-zero diagrams are 




1" 



4+ 




T)—^ [L] 3+ 4+ {l 

P41 ^ 



\ 




The relevant lower-point amplitudes with a (p can be found in [20] 
A\(P,l-,2-) = ^A'{cp,l~,2-) 

R{cP, 1-, 2-, 3+) = 1", 2", 3+) 

AU6 I- 2+ 3+) - ^ f 1 - ^ + ^^'l 
A [<p,l ,2 ,3 )- 43^2^1 n) (23)2 

-^A0(</)t,l-,2+,3+) 

while the last was calculated in [46] 

AUl~ 2+ 3+ 4+) = ^ f 1 - ^ + ^] 
^ '-^ ) 48vr2V iVy[12](23)(34)[41]' 

This gives the direct recursive terms 

DR23 = A3 (2- , 3+ , -P23 ) T^r ^4 4+ , T- , P+ ) 

^ cr / n/ , [34]gi23[13](4Pi23] 
3 V Nj [12][23](4P23l]2 

-2cM°((/.t,l-,2-,3+,4+) (8.48) 

L'i?234 = i?4(2-,3+,4+,-p2^J-^^3(0, 1^,^234) 

-^234 

= _!V , M j34]s|234M!iMklL f84m 

3\ iV^iV;(34)2s234(2P34l](4P23l]2 ^' ^ 



.46) 
.47) 
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DR^i = A3 (4+ , 1 - , - p+ ) J^i?4 (0, 2- , 3+ , ) 

-Ml 



.50) 



L>i?34i = i24 (3+ , 4+ , 1- , ) -^^3 (0, 2- , P3-41 ) 

-^341 



iV + ivj(34).S34i[41]2(2P34l]' ^^'^^^ 

To obtain i?, we now have to add the DR terms and CR and subtract the O 
terms. Some lengthy but trivial manipulations give 

R = 2cr^A^A, 1", 2" , 3^, 4^) + Ci^ + |M (l - ^ + ^ 
(23)(1P243]2 (41)(3Pi23] ^ (14)(2Pi34]2 



(34)[43][32]s234 (34) [12] [32] (34)[43][41]s34i 

(32)(4Pi24] (12)(2Pi34] (12) (1^243] 

(34) [12] [41] 2[41]s34i 2[32]s234 
-(12)2 — ^ — ^ — ^ 



S34 S12 2S23 2s, 



41 



^.52) 



It is of special interest here that the first term is present regardless of super- 
symmetry, where we might have expected that an amplitude with supersymmet- 
ric particle content running in the loop would be fully cut-constructible. The 
catch is, of course, that cut-constructibility requires also the external particles 
to be supersymmetric, something which is not the case with (p. With those ar- 
guments it may seem slightly surprising (although welcome) that it drops out of 
^i(i/, l-,2-,3+,4+) when we use (|8.11|) . In A^iA. 1 ,2 , 3^, 4"*"), of course, the 
additional rational piece stays. 



8.2.3 Tests 

There are a series of tests which we can perform to check that the result is correct. 
The first regards only the cut-containing part, in particular the IR divergences 
in e. These are given by (c/. section [5.1.4p 

^1 = -££^0 V ( + 0{e\ (8.53) 

This can be checked readily by taking the IR divergent parts of the box and 
triangle functions from appendix |A] and putting them into ()8.33p . There are no 
contributions from the bubble- like and rational parts. 
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CoUinear Limits 

The second test is that of cohinear Hmits. The full checking of all collinear limits 
is a rather tedious exercise and will not be repeated here. It is contained in 
the original article in detail. Since the cut-containing parts have been fixed by 
unitarity, we expect the right collinear limits automatically; for the rational parts, 
the use of recursion amounts exactly to enforcing the correct collinear limits for 
gluons 4 and 1, and for gluons 2 and 3. Checking these limits provide primarily 
a check on the algebra. 

From the perspective of on-shell recursion, the important collinear limits are 
the rational ones which do not participate in the direct recursive terms. If these 
turn out wrong, it is a sign that the use of recursion was not justified because 
R{z) — CR{z) did not vanish as z ^ cxd. In the case at hand, the important 
limits in this sense is when 1 and 2 become collinear and when 3 and 4 become 
collinear. 

For these checks, we can use the rational parts of the splitting functions for 
a scalar running in the loop (j5.39p along with the tree splitting functions (|2.52p . 
These predict that 



as 3 and 4 go collinear. We can prove the first by first extracting the part of 
(|8.52|) that has the correct singularity. 





as 1 and 2 go collinear, and 





ni n£\ r (41)(3Pi2 3] (32)(4Pi24] 
N n)[ [12] [32] [12] [41] 



-1 



2ctA 



(0t,p-,3 






-1 
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as wanted. For the colhnear hmit between 3 and 4 things are a bit more messy. 
Because R{(t),l- ,2- , P+) ~ A°{ct>,l' , P+), the first hne of (1835]) comes from 
the first term of (18.520 in the same way as for 1 and 2. The remaining terms with 
the colhnear singularity are 



cr[34] / _ ^ , ^ 
3(34) V N N 



(23)(1P243]^ 



+ 



(14)(2P3i4]^ 



(41)(3Pi23] 
(34) [12] [32] 



cr[34] 
3(34) 



n 



n. 



1 - ^ + — 

N N 



(34)[43][32]s234 (34)[43][41]s34i 
(32)(4Pi24] (12)^ - 
(34) [12] [41] (34)[43]_ 

(23) (12) [23] ((12) [23] +2(14) [43] 



(34)[43][32]s234 



(14)(21)[14]((21)[14] +2(23)[34] 



(12)^ 



(34)[43][41]s34i 
y^z{l-z) f{Pl){PPi2P] 



cr[34] 
3(34) 



cr[34] 
3(34) 



(34) 



N N 



[12] [P2 
(12)2" 



+ 



(34) [43] 

(P2)(PPi2P] 



S23 



+ 



[12] [PI 

S41 



S34 ^ S234 ■S34I 



1+2 



(12) (23) (41) 
(34) 



1 



+ 



1 



S234 S341 
N N 



+2 



^/z{l-z) {PP12P? 

(34) [12][2P][P1] 
- (12)3 [24] [13] ^ (i2)2[2l](41)(23) 

[34]S234S341 (34)S234S341 

(12)(23)(41)((3Pi2 3] + (4Pi24]) ./^l 



c^^/zil 



(34)S234S341 



N N 



(34) 



(12) 



[12][2P][P1] 
[34]m^ 



(34)(2P)(P1) (34)2[12][2P][P1] 



n 



n. 



1 - ^ + ^ 



N 
,2 



N 
P+) 



(34) 



+ 



[34]A^{^,l-,2-,p-) 
(34)2 



^.57) 



as wanted. Notice how we had to rewrite a contribution which superficially went 
as s^^ to one with either a (34) ^ or a [34] -1 diver gence. As noted above, the 
correctness of these two colhnear limits are an indicative necessary condition for 
the recursion relation to hold. 



Soft Higgs Limit 

The very last test we will make is the limit as the Higgs momentum goes soft. In 
fact, it turns out to fail, but only because the naive expectation for the soft limit 
is wrong. The naive expectation comes from noting that as the Higgs momentum 
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goes to zero, the H field in the effective correction CHF^^F^^'^ becomes a constant, 
and the term becomes proportional to the Higgsless gluon term. This tells us that 
the first order term in C should be obtainable in this limit as 

d 

AHkn 0,n gluons) = Cg—AHn gluons) 

dg 

= C{n- 2 + 21) A\n gluons), (8.58) 

where g is the gauge coupling and / is the number of loops. From looking at the 
form of the (/)-MHV amplitude and the (/)^-googly-MHV amplitude, it is reasonable 
to guess that the corresponding limits for (f) and (j)^ are 

A\k^ 0, n.g-,n+g-^) = (n_ - 1 + l)A\n_g- ,n+g+), 
A\k^t ^ 0, n.g~,n+g+) = (n+ - 1 + l)A\n^g- ,n+g+). (8.59) 

In other words, we should be trying to prove that 

^i(</.,l-,2-,3+,4+) ^ 2^1(1-, 2-, 3+, 4+) (8.60) 

as k(j, 0. ^^(1^,2^,3^,4+) in the FDH renormalization scheme has the cut- 
containing and rational parts [46] 



C(l-,2-,3+,4+) = crA"{l-,2-,3 



2Fl""(l,2,3,4) 



--^-^V2(-23) 
3 3N 3N J ^ ' 

i?(l-,2-,3+,4+) = - ^ + ^)^°(l-,2-,3+,4+). (8.61) 

The cut-containing part can be checked quite easily by taking the expression 
(I8.33P and converting the functions back into F2 bubble functions and then 
treating each term in turn. In short, only the one- mass box (which becomes 
zero-mass) and -F2(s23) survive the limit, and for the bubble functions it is only 
the last line of (j8.33p that survives. 

The limit of the rational part is in principle done in the same way. The term 
proportional to ^'^(A, 1~, 2~, 3+, 4+) vanishes in the limit, and there are several 
terms which must be considered together if the limit has to be taken consistently. 
In the end we find that 

i?(<^,l-,2^3+,4+)^ -^(^1-^ + ^)^0(1-, 2-, 3+,4+), (8.62) 

as /C(^ 0. The corresponding 0^ rational part has the soft limit 
i?(</.t,l~,2-,3+,4+) = i?(0,l+,2+,3-,4-)* 

1 



487r2 



N N J ^ 



l--rf + ^ U"(l+,2+,3-,4-)* 
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such that the Higgs soft hmit is 



i?(iJ,l-,2-,3+,4+)^-^( 



which is a factor —1/3 away from tlie expectation. Tliis sort of violation was 
anticipated aheady in [20] where it was pointed out that the Higgs momentum can 
act as a kind of effective infrared regulator and that this would cause an exchange- 
of-limits problem with the dimensional regularization. Thus, we conclude that 
our results are indeed correct. 
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Appendix A 

Integral Functions 



The box, triangle and bubble integrals used in this thesis can be obtained from 
[34], though with slightly different conventions. The boxes are given by 



f d'^-'^'L 1 



—i/j, 



2e 



(A.l) 



where the usual ie prescribtion is understood. The integrals are classified ac- 
cording to which of the momenta are massive and which are masslcss. 1-mass 
integrals have one massive momentum, 4-mass integrals have four. There are two 
different 2-mass configurations; if the massive momenta are opposite it is called 
the 2-mass-easy and if they are adjacent it is called the 2-mass-hard. 

We give the integrals in terms of the box functions i<4(. . .) which are dimen- 
sionless because they are multiplied by the Gram determinant Q of that integral. 
We also take out a factor rr coming from dimensional reduction. 



cr 



where 



cr 



1 r(l + e)r2(l-e) 



(47r)2-f r(l - 2e) IGtt^ 
We also define the kinematic invariants 

s = {Ki + K2f, t = iK2 + Ksf. 
The box functions through 3-mass are. 



+ 0{e). 



(A.2) 



(A.3) 



(A.4) 



1 



,2 \ e 



+ 



2 \ e 



+ log2 



-t 



(A.5) 



Fr{ki,k2,k3,K4) 
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-Ki 



-Li2 1 



+ Lio 1 



1 



2 \ e 



/i 



+ 



2 \ e 



+Li2 1 



Lio 1 



st 



2 \ e 



2 \ e 



+ 1^1 - 



-^1 



1 



1 



2 \ e 



— s 



+ 



2 \ e 



-^1 



-^1 



1 



KiKl) 



+ 



KlKl 



+ Li2f 1 



The Gram determinants are given by 



gOm ^ gim ^ g:.ma ^ 

Triangles are defined by the integral 

2 f d^-^^L 



2mh 



(A.6) 



(A.7) 



(A.8) 



(A.IO) 
(A.ll) 



J W: 



(A.12) 



again with ie understood. They have the same classification in terms of massive 
corners as the boxes. As in the box case, we can take factors out of the integral 
and define triangle functions 



(A.13) 
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The triangle functions are (omitting the 3-mass) 



-Kl 



where 



glm 
Q2m 



-Kl + Kl 



The bubble integral is given by 



hiK) = -ifi 



2t 



1 



-Kl 



cr 



(2^)4-2^ L2(L + Ky 

2 



e(l-2e) \-K^ 
= -r[- + log^-^ 
= crF2{K^). 



+ 2 



+ Oie) 



(A.14) 
(A.15) 



(A.16) 



(A.17) 
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Appendix B 



Solutions of Quadruple Cut 
Constraints 



This appendix lists the solutions to the cut constraints in a quadruple cut of a 
one-loop amplitude. For each number of massive corners there are two solutions 
which are distinguished by which massless corners are MHV and which are googly- 
MHV. The two solutions will always be related by flipping (|-) |-]) but they 
are all included for clarity. The four-mass cut is not used in this thesis and the 
reader is referred to the original article for those solutions [59] . It should also be 
noted that all the solutions here follow from the three-mass cases (which follow 
from four-mass) by taking external momenta lightlike, but again, all cases are 
included to benefit maximally from the masslessness of some corners. 

Our conventions are to denote the internal on-shell momenta h through I4 
starting from the bottom of the diagram and moving clockwise around the loop. 
These momenta point in the clockwise direction. The corners are denoted by 
A through D starting in the lower left corner. The outgoing momenta at those 
corners are denoted by the corresponding letter in uppercase when it is massive 
and in lowercase when it is massless. Massless corners are denoted by © when 
they are MHV and by © when they are googly-MHV. 




\d){dCBA\ 

{dACd) 
\BCd){dA\ 

{dACd) 
\BAd){dC\ 

{dACd) 
\d){dABC\ 

{dACd) 



(B.l) 
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B 



A 
b 



A 
b 



A 



B 




A 




h 
h 
h 
h 

h 
h 
h 
k 

h 
h 
h 
U 

h 
h 

h 
k 



ABCd][d\ 



[dCAd] 
Ad][dCB\ 



[dCAd] 
Cd][dAB\ 



[dCAd] 
CBAd][d\ 



[dCAd] 



d){bA\ 



(bd) 
b){dA\ 



(bd) 
b){dC\ 



(db) 
d){bC\ 



{db) 



Ab][d\ 



[db] 
Ad][b\ 



[db] 
Cd][b\ 



[bd] 
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Appendix C 



Passarino— Veltman Reduction 
of Cuts 



This appendix deals with the reduction of cut integrals which have loop momen- 
tum dependence in their numerators. Those needed for this thesis are: 

Io = Jdn, I^ = Jl''dn, I^" = J n^dn, (C.l) 

ji^= j dii r-=[-^df, (C2) 

J {mlm] ' ''■ J {mlm] ' ^ J {mlm] ' 



where 



dfi = 6{1^)6{{1 - PY)dH. (C.4) 

The goal is to express all of these integrals in terms of m'^, P^, g^'^ , Iq, and Jq. 
The method can be illustrated for by noting Lorentz invariance forces it to 

take on the structure 

jM ^ pMj' + ^Mj"^ (C.5) 
where /' and /" are scalar integrals. These can be determined from the equations 

(mJim] = [ ^j^^^dn = Io = imPmU' (C.6) 
J {mlm\ 

and 

/op . / 

— — ^d/x = P^Jo = 2P^I' + (mPmU". (C.7) 
{mlm\ 

Similarly, II^'^ can be computed by noting that it must take on the form 

I^'' = gt^'^l' + pMp^/" (C.8) 
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and using the contractions 

I^^ = J fdn = = 4l' + P^I", (C.9) 

and 

2pi/jM'^ = J{2P- l)l''dn = P^I'^ = 2P''I' + IP^Pf^l", (C.IO) 

which relate /' and /" to This procedure will allow us to construct the 
required integrals recursively. Rather than deriving the results in detail, we will 
just note that they are 

if = (C.ll) 
= -^h + \p^P''h (C.12) 

^ {mPm] " ~'~ \ {mPmY ^ ~^ {mPm] ^ ) ^ ^ 

ju,,^ _ g^'^P'^ J P^P" J {m^'P'' + P''m'')P'^ ^ 
^ A{mPm] ^ 2{mPm] ^ 2{m,Pm]'^ ° 

\ {mPm]^ ^ (mPm]^ '^y ' ^ ^ 

12(mPmJ 3(mPmJ 

^pUpU^K _^ ^fipUpK pH^KpU^p2 ^ 

6(mPm]^ ° 
(jnt^jni^P'^ + P'^m^m" + mt'P''m'')P'^ 
+ 3(;;^PH^ ^° ^^-^^^ 

(g^"-/»"-+.g'-";»/'+g''/'/»'OP% 

12(mPm]2 ^^-^^^ 
11 P6 , P6 



- y ^;^/o + ^^^,Jo ) . (C.19) 
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